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ABSTRACT 




The calculation model performs a rapid preliminary pitchline 
optimization of axial gas turbine annular flowpath geometry, as well as 
an initial estimate of blade profile shapes, given only a minimum of 
thermodynamic cycle requirements. No geometric parameters need be 
specified. The following preliminary design data are detenvned: 

(1) the optimum flowpath geometry, within mechanical stress limits; 

(2) initial estimates of cascade blade shapes; 

(3) predictions of expected turbine performance. 

The model uses an inverse calculation technique whereby blade profiles 
are generated by designing channels to yield a specified velocity distribu- 
tion on the two walls. Velocity distributions are then used to calculate 
the cascade loss parameters. Calculated blade shapes are used primarily 
to determine whether the assumed velocity loadings are physically 
realistic. Model verification is accomplished by comparison of predicted 
turbine geometry and performance with an array of seven NASA single-stage 
axial gas turbine configurations, ranging in size from 0.6 kg/s to 
63.8 kg/s mass flow and in specific work output from 153 J/g to 558 J/g 
at design (hot) conditions; stage loading factor ranges from 1.15 to 
4.66. 
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NOMENCLATURE 


critical sound speed, m/s 
lift coefficient 
diffusion coefficient 

2 

gravity constant, m/s 

channel or blade height, m 

rotor tip clearance gap. m 

velocity loading parameter defined by Eq. (23) 

distance along axial direction, m 

mean chord length in axial direction, m 

2 

static pressure, n/m 

difference between blade pressure-surface static pressure and 

2 

suction-surface static pressure, n/m 

2 

total or stagnation pressure, n/m 

radius coordinate, measured from turbine axis of rotation, m 

acceleration parameter, defined by Eq. (21) 

blade pitch or spacing, m 

blade thickness, m 

rotor speed, m/s 

velocity measured relative to non-rotating reference frame, m/s 

blade loading velocity defined by Eqs. (24a-b) 

velocity measured relative to rotating reference frame, m/s 

mass flow rat;., kg/s 

total or stagnation temperature, °K 

loss coefficient (pressure) 

number of blades 


a = flow angle measured with respect to absolute (non-rotating) 
velocities 

3 = flow angle measured with respect to relative (rotating) velociti 

y = ratio of gas specific heats 

6 “ pressure ratio referenced to standard conditions 

e « channel thickness parameter, defined by Eq. (18) 

n = adiabatic efficiency 

<)> = angular coordinate 

3 

p = density, kg/m 

fil - distance along axial direction normalized by chord length M 
6 = temperature ratio referenced to standard conditions 

a = cascade solidity 

Subscripts 
c = centerline 
eg = clearance gap 
e = cascade exit 

i = cascade inlet 

m = maximum 

M = arbitrary location along axial direction 

p = pressure surface 

PA = profile loss 

s = suction surface 

sf = secondary flow loss 

std = standard conditions 

tot = total 


viii 


> 

^ u = tangential 

x - axial 

Superscripts 

* 1 s absolute total or stagnation condition 

" - relative total or stagnation condition 

* = parameter normalized with respect to conditions at cascade inlet 

> 

) 

> 
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INTRODUCTION 


Axial gas turbine design can be described as a search for an 
optimum set of geometric parameters to satisfy certain general specifica- 
tions such as work output and mass throughflow. Traditionally, the 
search procedure begins with computations which provide a quick estimate 
of optimum overall configuration for a given cycle data point. Such 
computations usually rely upon loss correlations that mainly depend on 
cascade inlet and exit conditions, without regard to specific cascade 
geometry, t.'se of such correlation models implicitly assumes that 
detailed blada geometry does not significantly influence flow losses, 
an assumption which is open to question. 

The relationships between blade shape, loading, and flow losses 
are established in more detail later in the design procedure, a process 
which requires detailed knowledge of individual blade shapes. Since the 
overall design procedure is an iteration loop involving geometry and 
aerodynamic losses, however, these relationships should be established 
as early in the preliminary design procedure as possible. This should 
result in a more realistic final design in less overall time. 

The present optimization technique seeks to provide an initial 
axial turbine stage design procedure which directly links overall stage 
performance and flowpath geometry with internally generated cascade 
loadings and blade shapes. Having such detailed geometry generated 
internally has a threefold advantage: 


(i) excessive data preparation times are eliminated; 

(ii) the loss model is provided flexibility in determining blade 
shapes (loadings) which optimize boundary layer (Reynolds 
number) effects; 

(iii) relations between such factors as blade chord, blade camber- 
line, blade stagger, etc., and the specified velocity 
triangles are automatically accounted for without additional 
correlations. 

The present design procedure is "preliminary" in that calculations 
are performed only for an average streamsheet surface (pitchline) location 
within the turbine stage. Pitchline radius may vary (linearly) within a 
blade row; however, radial velocity components are neglected and such 
variation should be small. The most basic assumption of the analysis is 
that velocity loadings on the blading are known. These loadings are used 
to calculate blade profile shapes and profile (friction) losses. The 
blade shapes are, in turn, examined to determine whether the assumed 
loading is physically realistic. 

CALCULATION OF CHANNEL (BLADE) SHAPES 

Cascade channel shapes (blade profiles) are generated using a 
simplified technique to provide solutions to the so-called indirect 
problem: that of the design of a channel to yield a specified velocity 

distribution on the two walls. Adjacent walls then become the suction 
and pressure surfaces, respectively, of the individual blades. The 
technique is similar to that utilized by J. D. Stanitz^. 
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For given thermodynamic cycle requirements, wheel speed, and 
adiabatic efficiency, stage velocity triangles can easily be obtained 
along some representative pitch-line radius. If the channel velocity 
distribution is also known, the channel (blade) shape is generated as 
follows: 


Continuity 


w . cos a, (2,rR i h, ) (p'^ r > std 

p ®cr . * y i 


( 1 ) 


Moment of Momentum 
W 

gz^; (RV u )dm = (4 P) hRdl " (2) 

Integration of Eq. (2) along the entire axial chord length of the blade 
results in the overall change in tangential (swirl) velocity across the 
cascade: 

1 

(RV u ) e - (RV U ) 1 = £ h i R 1 .MS,P std { h*R* ff dfl (3) 

r 0 1 

Note that "axial" chord is actually the vector sum of axial and radial 
displacements through the cascade, since the pitch-line radius may vary. 
Making use of the continuity equation (1), and the cascade solidity, a, 
defined as the ratio of axial chord to blade spacing, Eq. (3) can be 
rearranged as 
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If the integration of Eq. (2) is carried out only to some 
arbitrary point fl] within the blade row, the result is 
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from which the local absolute tangential velocity at that point, . 
<V a ir> 01 can be obtained. A more useful relation is the tangential 
velocity distribution relative to a blade row moving at some wheel speed 
(U/a^ r )^. defined by 
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(7) 


The term a' is the local critical velocity defined with respect to 
cr 

absolute total temperature. Now, since critical velocity can also be 
defined with respect to relative total temperature, Eq. (6) can be 
rewritten as 
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or upon combining Eqs. (5), (6), (7), and (8), 



Eq. (9) specifies the relative tangential velocity component at any 
location within the cascade. 

GENERATION OF BLADE CAMBERLINE DISTRIBUTION 

At this point in the computation it is possible to analytically 
define a mean camberline distribution for the cascade. Let W/a£ r and 
W/a" represent the relative critical velocity ratios along the 

Ci |S 

pressure and suction surfaces of the blade, respectively; note that both 
quantities are known since the velocity loading is presumed to be known. 
It is desired to calculate the average flow direction, 8^, for any 
location fil within the blade passage. In order to do this, an additional 
assumption is required. It might be assumed that the channel mean 
velocity is the arithmetic mean of the suction surface and the pressure 
surface velocities. 
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On the other hand, one might assume, say, a linear variation of axial 

critical velocity ratio, W /a" , through the blade passage, in which 

cr 

case 




= tan 
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and (W/a£ r )^ is approximated by 


(W/aJr)^ • 
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( 13 ) 


Even more flexibility can be provided by modifying the linear variation 
by a variable increment; the increment is zero at the cascade inlet and 
exit, and some specified maximum value at a given location within the 
cascade. The increment must then satisfy four boundary conditions and 
is defined by a third order polynomial. Experience has shown that either 
the linear or modified linear axial velocity distributions provide the 
most realistic blade shapes. 

If the mean camberline is described by the cylindrical coordinates 
R, m, and tj>, then the camberline is generated by the expression 
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3m f^) " R ff + * 3i = tan 8 + * 3i (,4) 

Equation (14) must be integrated point-by-point through the cascade. 

GENERATION OF LOCAL BLADE THICKNESS 

For any location ft along the cascade passage the continuity 
equation can be written as 

W f 

= (pW cos 6) (h x channel width) 

Define a "channel width parameter," e, by 

- channel width _ channel width 
e blade pitch 2 ttR 

z 

Equation (15) then becomes 
Wf = (pW cos 3 ) ( 2TrehR) 

for a given position /fl. Combination of Eqs. (1) and (17) yields, after 
rearrangement. 
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Now local blade thickness, t^, is the difference between local blade 
pitch and local channel width, so that Eq. (16) can be written as 
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Local blade thickness is then distributed equally about the camberline 
distribution defined by Eq. (14). Specification of the blade profile 
is thus complete. A representative blade section generated by this 
procedure is shown in Figure 1. 



Figure 1: Representative Blade Section Generated by Inverse Method 

CALCULATION OF BLADE LOADINGS 

Of fundamental importance to the present optimization procedure 
is the selection of proper cascade blade loadings. These loadings are 
used to 1) calculate channel shapes (blade profiles), 2) meet required 
cascade solidity specifications, and 3) calculate profile (friction) 
losses. The loading model should be simple, yet flexible enough to 
acconmiodate a wide range of channel centerline accelerations, blade 
surface diffusion rates, and cascade loading requirements. 

The loading model proceeds from the definition of a centerline 
velocity distribution, V ’ , as 

V 


» 
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1 / 0 ! 
p-y Esin tt(/9 - I) - 1] + 1 


(20) 


The parameters V Q and V represent (dimensionless) critical velocity 
ratios measured relative to the blade row, with V measured at the exit 
of the row. The parameter ft is indicative of the channel centerline 
flow acceleration and is defined as 


ft = 1 


V. 

i 

TT 

v e 


(21) 


where I/.. is the critical velocity ratio relative to the blade row, 

measured at the inlet to the row. It should be noted that 1/ is a 

c 

fictitious parameter used only in the definition of blade loading and 
is not necessarily related to channel mean velocity, Eq. (13). 

Suction Surface Velocity 

It is assumed that the behavior of the velocity distribution on 
the blade suction surface is of primary importance to the overall 
optimization process. Although velocity diffusion occurs on both blade 
surfaces, the adverse effects of pressure-surface diffusion are to some 
degree obviated by the subsequent downstream re-acceleration of the flow 
along that surface; such re-acceleration does not usually occur on the 
suction side of the profile. 

A fundamental characteristic of the suction-surface velocity 

) 

distribution is deceleration of the boundary layer flow near the trailing 
edge of the blade row. Of interest is both the amount of such decelera- 
tion and the location along the suction surface at which such deceleration 
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begins. If velocity "spike" effects near the leading edge of the blade 
are ignored, velocity diffusion on the suction surface becomes a simple 
question of where maximum velocity occurs and what its value is. 

Following standard practice, a "diffusion coefficient" for the 
suction surface velocity is defined as 

1 / 

D s = 1 ' JT (22) 

where V m represents the maximum velocity occurring anywhere on the suction 
surface. Figure 2 illustrates each of the pertinent velocities employed 
in the model. As in Reference 2, the suction surface velocity is 
described by a simple piecewise parabolic curve fit; the point of 
maximum velocity, is also the point where the two curve portions are 
matched. It has been suggested that the location of is dependent on 
the value of the channel centerline flow acceleration ft as 

- 0.5 + 0.3[1 - (1 - J?) 2 ] (23) 



Figure 2: Blade Velocity Loading Model 


11 


Barring user input to the contrary, £ is calculated from Eq. (23). 
Should the user so desire, however, any value for £ m can be specified. 

The suction surface velocity distribution is defined by the 
following relations: 
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It should be noted that it is, at least mathematically, possible to 
have negative values for the diffusion coefficient; in this case, 

1/ < V n and there is no real "maximum" velocity (other than 1/ ). For 
this situation, the piecewise curve fit degenerates to a single 
second-order curve of the form 
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Pressure Surface Velocity 

The pressure surface velocity distribution is obtained by assuming 
that the entire blade loading is distributed equally about the centerline 
velocity distribution 1/.. The pressure surface distribution is then 
given by 


where is given by either Eq. (24) or Eq. (25). Though there is no 
mathematical constraint, negative (reverse direction) velocities are 
not desirable, and are not allowed by the model. 

OPTIMIZATION OF THE DIFFUSION COEFFICIENT 

The blade loading diagram, suction surface diffusion coefficient, 

3 

and cascade solidity can be tied together through the Zweifel loading 
coefficient. The Zweifel coefficient relates cascade solidity to flow 
angles at the cascade inlet and outlet, and, according to reference 3, 
tends to assume a certain narrow range of values for optimum cascade 
performance. Normal design practice utilizes Zweifel coefficients in 
the range of, say, 0.7 to 1.0. Thus, for specified flow angles and a 
given Zweifel coefficient within the "optimum" range, an optimum cascade 
solidity can be calculated. Now, solidity is related to the individual 
blade velocity loading through Eq. (4). Furthermore, for known cascade 
inlet and outlet conditions, the assumed velocity loading model possesses 
only one independent variable: the suction surface diffusion coefficient. 

Thus, 

(1) specified thermodynamic cycle requirements establish cascade 
inlet and outlet conditions; 

(2) these conditions, together with an input Zweifel coefficient 
will determine an optimum solidity; 

(3) the present model then adjusts the suction surface diffusion 
coefficient so that the velocity loading, when integrated 


i3 

in Eq. (4), yields a calculated solidity equal to the 
optimum solidity. 

The final velocity loading is then used to calculate both blade shapes 
and blade profile losses. 

CALCULATION OF CASCADE PERFORMANCE LOSS 

Three aero thermodynamic loss mechanisms are accounted for in the 
present analysis: profile loss, secondary flow loss, and rotor tip 

clearance loss. 

Profile loss is defined here as a combination of frictional 

effects resulting from the flow of a viscous fluid over a solid surface 

and the subsequent downstream mixing of the suction surface and the 

pressure surface boundary layers. Both losses are accounted for through 

4 

use of the Stewart mixing loss theory , which defines them in terms of 
overall displacement and momentum thicknesses. Secondary flow losses 
are due to the annulus wall boundary layers and their interaction with 
blade rows. Dunham's review paper presents an excellent analysis of 
the phenomenon. Dunham states that two separate effects must be 
accounted for if losses are to be properly estimated: 

(1) a vortex core loss, arising from fluid, originally in the 
upstream wall boundary layer, being subsequently shed from 
the trailing edge of the cascade; 

(2) a downstream wail boundary layer loss, wherein fluid origi- 
nally in the mainstream is entrapped in a boundary layer 
developing on the annulus walls within the cascade. 


Rotor clearance, or tip leakage, losses effect turbine performance 
in two ways: 

(1) leakage diminishes overall work extraction from the flowing 
gas, since some fluid at the tip is not turned, and 

(2) leakage produces an underturning of the flow felt in regions 
other than the tip region, which further diminishes work 
extraction. 


Profile Loss 

When calculating profile boundary layer losses it is possible 
to utilize models of almost any complexity imaginable. Of primary 
importance, and of least certainty, is the location along a given blade 
surface of the point of transition of the boundary layer from laminar to 
turbulent flow, as well as the point of separation of the boundary layer 
from the wall. Since the present analysis is' of a preliminary design 
(pitch-line) nature, exceedingly complex loss models are not appropriate. 

In its simplest sense, the location of the boundary layer transi- 
tion point is a cascade Reynolds number effect. As Reynolds number 
(based on cascade exit velocity and blade mean camberline length) is 
increased, the transition point tends to move upstream towards the leading 
edge of the blading. For the laminar portion of the boundary layer, the 
present analysis calculates loss parameters using the Truckenbrodt 
approximation for boundary layer growth on a wall with pressure gradient, 
as described, for instance, in Reference 6. For the turbulent portion, 
loss parameters are calculated using a simple formulation described in 
Reference 7. A rough estimate of transition is made using a "critical 
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Reynolds number" parameter, based on boundary layer displacement thickness, 
again described in Reference 6. Even though this procedure might be 
considered "too simple," it provides a way to account for Reynolds 
number effects on overall cascade loss, and is consistent with the 
fact that velocity loadings are already assumed to be known. As with 
any set of assumptions, the final justification lies with how well the 
model predicts reality. Pertinent relationships are described in 
Appendix B. 


Secondary Flow Loss 

The correlation for secondary flow loss as suggested in the 
conclusion of Dunham's paper requires a knowledge of the boundary layer 
displacement thickness on the endwalls upstream of the cascade, a parameter 
which is not readily available, at least with any degree of confidence 
(especially for the rotor). Actually, the upstream boundary layer, or 
vortex core, loss adds to a second loss component, called the downstream 
loss. The present analysis assumes that the sum of these two loss 
components can be represented as a constant value. The specific 
expression used for secondary flow loss then becomes 


Y 


sf 


0.0138 


M cos a e 
h cos a. 



cos 2 a e 

3 

cos a_ 
m 


(27) 


where is a pressure loss coefficient, is a cascade lift coefficient 
(defined only in terms of inlet and outlet relative flow angles), and 
the value 0.0138 represents the sum of the two loss components discussed 
above (see Figure 9 of Reference 5). The remaining parameters are 
described in the Nomenclature. 


Q 


Tip Leakage Loss 


Q 

Though recent work such as that of Lakshminarayana provides 
insight into the physics of the tip leakage loss mechanism, quantitative 
predictions of performance loss still require some correlation with 
experimental data. In this regard, the correlation reported in Reference 
7, taken from data originally reported in Reference 9, shows promise. 

This correlation indicates that, for given blade reaction, the tip- 
clearance loss varies (approximately) linearly with clearance gap: 
furthermore, the loss increases for increasing blade reaction. The 
latter is intuitively correct, since high reactions (large pressure 
differences) cause more high-kinetic-energy flow to leak through the 
clearance gap. 

The various curves given in Reference 7 can all be approximated 
by the single expression 


HfiSL = 
n 


(2.755 ft 2 + 0.108 /? + 1.72) 


where n cg 3 efficiency with t-'lea ranee 

n = efficiency without clearance 
h = blade height 
Ah = clearance gap 

and ft is blade reaction at the tip , defined by 


ft 




2 2 
- W. + V • 
l 1 


(28) 


( 29 ) 


where W, V represent relative and absolute velocities, respecti vely. 
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AERODYNAMIC OPTIMIZATION OF BLADE CHORD 

Both blade profile (friction) losses and cascade secondary flow 
losses are dependent on blade chord. Normally, blade chord is determined 
from manufacturing trade-offs or prior design experience. For instance, 
chords below a certain minimum length cannot be made economically; 
chords which are too large increase overall engine length (weight), 
etc. While chord lengths can be input to the present model, a 
theoretically optimum chord can be calculated simply from a considera- 
tion of the profile and secondary flow loss mechanisms. 

4 

From the Stewert analysis a cascade absolute total pressure 
loss Pg/P^ ca*i be determined. A "profile loss coefficient," Y , can 
be defined as 

i - p;/p; 

V = p;/ppi - p 6 /p;j (30) 


where p e is the static pressure at the cascade exit. A "total loss 
coefficient" Y t t can be formed by 


tot 


■ V + v 


sf 


(31) 


where Y gf is defined by Eq. (27). Normally, profile losses are reported 
as a function of cascade Reynolds number, defined on the basis Of cascade 
exit velocity and blade chord length. In the present model, cascade 
Reynolds number is varied systematically by varying blade chord while 
holding all other parameters constant. This is possible since the 
velocity loading is defined with respect to a dimensionless length 



parameter (fl). Figure 3 illustrates how the different loss coefficients 
vary with changing rotor chord (mean camberline length). Note that, 
generally speaking, decreases with increasing chord while Y^ 
increases (blade height is held constant during the procedure). 

According to this simplified model there is an optimum Reynolds number, 
or optimum chord, which minimizes total aerodynamic loss (excluding tip 
leakage). The optimum value will, of course, change as cascade inlet 
and outlet conditions change. 



12 3 4 

Reynolds No. ( x 10 5 ) 


Figure 3: Typical Variation of Cascade Loss Coefficients with Reynolds 
Number 

OPTIMIZATION OF ANNULAR FLOWPATH GEOMETRY 

The model considers a range of annular (hub/tip) flowpath 
geometries, based on calculated mean flow properties along a pitchline 
radius. The minimum value for pitchline radius is defined by the condi- 
tion of zero acceleration, of the channel centerline velocity in the rotor. 
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The maximum value for pitchline radius may be set by limitations on 
blade tip speed or by aerodynamic considerations such as rotor limit 
loading. In any event, the loss model described previously is 
applied to a range of configurations within these limits in order to 
determine the one configuration having optimum aerodynamic performance, 
as characterized by adiabatic total-total or total-static efficiency. 

Hub and tip radii at the rotor exit section are obtained from user- 
specified axial Mach number and swirl conditions. A schematic flowchart 
of the entire calculation procedure is shown as Figure 4. 

ESTIMATION OF MECHANICAL STRESS 

Since the present analysis considers such a wide range of possible 
designs for a single set of cycle constraints, it is desirable to calcu- 
late preliminary estimates of mechanical stress conditions for each 
design. Rough estimates of both average tangential disk stress and 
blade root stress are made utilizing expressions found in textbooks 
on turbomachinery design (such as References 10 and 11). Blade root 
stresses are calculated using a taper correction factor of 2/3, which 
corresponds to a linearly tapered blade with a ratio of blade tip area 
to blade hub area of about 0.35. The disk half-area is considered to be 
a trapezoidal section; the disk rim load is considered to be the sum of 
1) the load due to the blades themselves, and 2) the load due to a blade 
attachment region, also taken to be a trapezoid. It should be noted that 
these stresses are directly dependent on both blade chord and blade 
solidity, both of which are optimized in the present analysis. It is 
impossible, therefore, to completely separate considerations of 
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Figure 4: Axial Turbine Design Optimization Procedure 
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aerothermodynamic optimization from considerations of mechanical 
integrity. It should also be noted that, while mechanical stress is 
calculated by the moilel , no design decisions are made on the basis of 
such calculations. 

VERIFICATION OF THE PERFORMANCE LOSS MODEL 

Verification of the performance loss model is accomplished by a 
process wherein test data, including test rig geometries normally calcu- 
lated internally, are input to determine whether actual measured test-rig 
performance can be predicted. Thermodynamic cycle requirements input 
represent air-equivalent (cold) conditions rather than design (hot) 
conditions. Quantities such as measured mass flow, rotor torque 
(specific work), and rotor exit swirl angle are utilized. 

Seven (7) NASA turbines have been selected to provide a means of 
comparison between aerodynamic performance (adiabatic efficiency) as 
predicted by the model and aerodynamic performance as measued during cold- 
air testing. These turbines are designated as 

Case 1: Low-cost civilian turbojet engine (Reference 12) 

Case 2: Research turbine for high temperature core engine 

application (Reference 13) 

Case 3: 12.766 centimeter tip diameter (solid blade configuration) 

turbine (Reference 14) 

Case 4: First stage of a 4£ stage fan-drive turbine (Reference 15) 

Case 5: Compressor drive turbine for a 75 KW automotive engine 

(Reference 16) 
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Case 6: Uncooled core turbine with high work output (Reference 

17) 

Case 7: Low-cost turbofan engine (first stage of a two-stage 

turbine) (Reference 18) 

These turbines range in size from 0.6 to 63.8 kg/s (1.3 to 140 lb/s) mass 
flow and in specific work output from 153 to 558 J/g (65 to 240 Btu/lb), 
all at design (hot) conditions; stage loading factor ranges from 1.15 
to 4.66. For comparison, relevant design parameters for each turbine 
configuration are given in Table 1. 


Table 1: Comparison of Relevant Turbine Parameters 
(Design Conditions) for NASA Test Cases 


Case 

Inlet 

Temp 

°K 

Inlet 

Pressure 

N/cm2 

Mass Flow 
Rate 
kg/s 

Specific Work 
Output 
J/g 

Work 

Factor 

1 

1089 

26.3 

3.28 

159.3 

1.15 

2 

2200 

386.1 

63.82 

287,3 

1.70 

3 

1478 

91.2 

0.95 

307.3 

1.67 

4A 

378 

24.3 

5,84 

25.7 

4.66 

5 

1325 

39.8 

0.60 

198.1 

2.10 

6 

2200 

386.1 

49. 41 

557.7 

1.94 

7 

978 

28.5 

2.99 

152.8 

1.72 


Equivalent design requirements— actual conditions not 
given. 


For each case, the following calculations are made: 

(1) With no geometry input, determine the optimum rotor exit 
mean radius by maximizing the turbine adiabatic total-total 
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efficiency. This calculation is performed with design (hot) 
conditions of thermodynamic parameters of flow, work, etc., 
as well as geometric design values such as rotor tip clearance, 
and is done merely to demonstrate the flexibility of the model. 
Thus, total-total efficiency as predicted by, say. Figure 5, 
need not compare on a one-to-one basis with that of Figures 
8 and 9. Further, n values shown in the geometric optimization 
portion of the present work differ slightly from those values 
given in Reference 19 due primarily to changes in the blade 
loading portion of the model. 

(2) With detailed test-rig geometry and flow conditions input, 
determine the predicted total-total efficiency as a function 
of 

(i) rotor solidity, a R 
(ii) rotor exit swirl angle 

In addition, for Case 1 only, a survey of optimum predicted flowpath 
geometry is made as a function of rotor Zweifel loading coefficient, 
and rotor exit swirl angle. 

Case 1^ 2 

Case 1 represents a turbine designed for a low-cost civilian 
turbojet engine application. The turbine configuration is single-stage, 
axial flow with a free-vortex whirl distribution. The aerodynamic design 
is conservative (low work factor) with moderate gas temperatures. Table 2 
presents the relevant design parameters and measures test-rig values 
input to the performance model . 


Table 2 : Case 1 Performance Model Input 


Parameter 

Design 

Equivalent 

Condition 

Test-Rig 

Condition 

Rotative Speed, RPM 

18075 

100% OF DESIGN 

Inlet Temperature, °K 

288.2 

310 

Inlet Pressure, N/cm* 

10.14 

10.8 

Specific Work, J/g 

43.23 

1.7% more 

Mass Flow, kg/s 
Rotor Exit Mean 

2.51 

9.5% LESS 

Radius, cm 
Rotor Exit Annulus 

10.06 

SAME 

Area, cm^ 

269.6 

SAME 

Rotor Clearance, cm 

0.028 

SAME 

Stator Solidity 

1.038 

SAME 

Rotor Solidity 

1.640 

SAME 

Stator Axial Chord, cm 

1.90 

SAME 

Rotor Axial Chord, cm 

1.76 

SAME 

Stator TE Thickness, cm 

0.101 

SAME 

Rotor TE Thickness, cm 

0.101 

SAME 

Rotor Exit Swirl, deg 

-3.8 

+14.0 


Figure 5 represents a preliminary optimization of annular flow- 
path geometry, wherein no geometrical constraints are input to the model. 
An optimum size for the turbine is chosen based on maximum total - 
total efficiency. The predicted size {rotor exit mean radius) compares 
well with the actual (design) size for the Case 1 configuration. 
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ROTOR EXIT 
MEAN RADIUS, cm 

Figure 5: Case 1 Preliminary Geometry Optimization Study 

Figures 6 and 7 demonstrate the predicted variation of optimum 
turbine size with rotor Zweifel coefficient (Figure 6) and rotor exit 
swirl (Figure 7), all other factors being constant. Figures 8 and 9 
are calculations wherein geometry is input to the model. Figure 8 is 
a parametric study of stage efficiency vs. rotor solidity (blade number); 
note that exit swirl is constrained to +14°, as opposed to the design 
value of -3.8°. At the design solidity, the model predicts an efficiency 
of 91.4, which compares well with the rig measured value of 91.0. 

Original design (target) efficiency was 88.0. Figure 9 represents a 
parametric study of predicted efficiency vs. rotor exit swirl, all other 
factors being constant. 
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Figure 6: Case 1 Geometry Optimization Study; Rotor Zweifel Coefficient 

Variation 



ROTOR EXIT 
MEAN RADIUS, cm 


Figure 7: Case 1 Geometry Optimization Study; Rotor Exit Swirl Variation 
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Figure 8: Case 1 Efficiency vs. Rotor Solidity; Comparison of Code Pre- 
diction, Original Design Va'lue, and Test-Rig Measurement 
(Geometry Input) 


13 

Case 2 IJ 

Case 2 is a half-scale model of a 50.8 cm (20 inch) turbine 
characterized by low aspect ratio, thick trailing edges, low solidity, 
and relatively large rotor tip clearance. Originally, the rotor 
blades had a constant section profile from hub to tip with no twist, 
resulting in relative ease of manufacture but a possible performance 
penalty. Subsequently, a free- vortex twist rotor was fabricated and 
tested. Due to the pitchline nature of the model used in the present 
code, these two conditions cannot be differentiated. Though the measured 
performance of both designs is reported (see Figures 11 and 12), only 
the untwisted rotor test-rig conditions are shown in Table 3. 


-15 - 10-5 0 5 10 15 20 25 


Figure 9 


ROTOR EXIT SWIRL, 
DEGREES 


Case 1 Efficiency vs. Rotor Exit Swirl for Fixed Stator and 
Rotor Solidities (Geometry Input) 





Table 3; Case 2 Performance Model Input 


Parameter 

Design 

Equivalent 

Condition 

Test-Rig 

Condition 

Rotative Speed, RPM 

12388 

100% OF DESIGN 

Inlet Temperature., °K 

288.2 

306 

Inlet Pressure, N/cm 2 

10.14 

17.24 

Specific Work, J/g 

39.57 

0.5% LESS 

Mass Flow, kg/s 

1.207 

1,5% LESS 

Rotor Exit Mean 


SAME 

Radius, cm 

11.75 

SAME 

Rotor Exit Annulus 


SAME 

Area, cm 2 

140.6 

SAME 

Rotor Clearance, cm 

0.043 

SAME 

Stator Solidity 

0.929 

SAME 

Rotor Solidity 

1.487 

SAME 

Stator Axial Chord, cm 

1,905 

SAME 

Rotor Axial Chord, cm 

1.715 

SAME 

Stator TE Thickness, cm 

0.089 

SAME 

Rotor TE Thickness, cm 

0,089 

SAME 

Rotor Exit Swirl, deg 

-17.8 

-11.5 


Figure 10 represents the preliminary optimization of annular 
flowpath geometry. As before, the predicted turbine size compares well 
with the actual size. 

Figure 11 is & fixed geometry parametric study of efficiency vs. 
rotor solidity, with all other factors held constant. Rotor exit swirl 
is fixed at the test-rig measured value of -11.5°. The measured 


TOTAL - TOTAL EFFICIENCY 
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ROTOR EXIT 
MEAN RADIUS, cm 


Figure 10: Case 2 Preliminary Geometry Optimization Study 
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ROTOR SOLIDITY, 

Figure 11: Case 2 Efficiency vs. Rotor Solidity; Comparison of Code 

Prediction, Original Design Value, and Test-Rig Measurement 
(Geometry Input; 

efficiency is shown as a band between 87.0 (untwisted rotor configuration) 
and 88.0 (twisted rotor configuration). Predicted performance is 86.6, 
compared to an original design value of 87.0. Recall that the present 
model cannot differentiate between the two rotor configurations. Figure 
12 presents a parametric study of efficiency vs. rotor exit swirl angle, 
with all other factors held constant. 

Case 3 1 4 

Case 3 represents an uncooled solid-blade version of a cooled 
turbine design in the 1 kg per second, 225-375 KW size class. Work factor 
and solidity are considered to be near optimum by conventional design 
standards. Both the stator and the rotor blading are untwisted and 
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ROTOR EXIT SWIRL, 
DEGREES 


Figure 12. Case 2 Efficiency vs. Rotor Exit Swirl for Fixed Stator 
and Rotor Solidities (Geometry Input) 

untapered. Typical of such small turbines, this design suffers from 
relatively large rotor tip clearance and secondary flow losses. Relevant 
test-rig conditions are given in Table 4. 

Figure 13 demonstrates that the predicted annular flowpath 
geometry agrees well with the actual geometry. Figure 14 represents 
a fixed geometry parametric study of efficiency vs. rotor solidity and 
indicates that the rotor solidity is indeed near optimum for the thermo- 
dynamic cycle requirements of case 3. The model predicts an efficiency 
of 34.1 for the design solidity compared to a measured efficiency of 
83.2. The original design efficiency was 85.0. Figure 15 presents a 
study of predicted efficiency vs. rotor exit swirl distribution. In the 
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Table 4: Case 3 Performance Model Input 


Parameter 

Design 

Equivalent 

Condition 

Test-Rig 

Condition 

Rotative Speed, RPM 

31460 

100% OF DESIGN 

Inlet Temperature, °K 

288,2 

300 

Inlet Pressure, N/cm* 

10,13 

8.27 

Specific Work, J/g 

62.1 

1.4% LESS 

Mass Flow, kg/s 

0.246 

6,1% LESS 

Rotor Exit Mean 



Radius, cm 

5.86 

SAME 

Rotor Exit Annulus 



Area, cm^ 

38,72 

SAME 

Rotor Clearance, cm 

0,025 

SAME 

Stator Solidity . 

1.098 

SAME 

Rotor Solidity 

1.551 

SAME 

Stator Axial Chord, cm 

0.721 

SAME 

Rotor Axial Chord, cm 

0.968 

SAME 

Stator TE Thickness, cm 

0.038 

SAME 

Rotor TE Thickness, cm 

0.050 

SAME 

Rotor Exit Swirl, deg 

-17.5 

-11.0 


range of swirl angles shown the total -total efficiency appears to be 
almost independent of swirl. The region denoted "model limit" represents 
choked flow conditions at the stator exit. 



ROTOR EXIT MEAN RADIUS, cm 


Figure 13: Case 3 Preliminary Geometry Optimization Study 

> 



ROTOR SOLIDITY, 

Figure 14: Case 3 Efficiency vs. Rotor Solidity; Comparison of Code 

Prediction, Original Design Value, and Test-Rig Measurement 
(Fixed Geometry) 
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ROTOR EXIT SWIRL, 
DEGREES 


Figure 15: Case 3 Efficiency vs. Rotor Exit Swirl for Fixed Stator and 

Rotor Solidities (Fixed Geometry) 


Case 4 represents the first stage of a 4£ stage turbine designed 
for high stage work factor. This turbine is characterized by shrouded 
rotors, high turning in both stator and rotor blade rows, and nearly 
symmetrical mean-radius velocity diagrams. As part of a development 
program, the first stage alone was fabricated and its performance deter- 
mined in cold air. Pertinent parameters are given in Table 5. 

Figure 16 presents the geometry optimization portion of the 
parametric study for Case 4. In this instance the predicted and actual 
sizes are not in good agreement; it should be noted, however, that this 
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> Table 5: Case 4 Performance Model Input 


1 

Parameter 

Design 

Equivalent 

Condition 

Test-Rig 

Condition 


Rotative Speed, RPM 

3098.7 

100% OF DESIGN 


Inlet Temperature, °K 

288.2 

378 

1 

Inlet Pressure, N/cm^ 

10.13 

24.3 


Specific Work, J/g 

25.65 

100% OF DESIGN 


Mass Flow, kg/s 
Rotor Exit Me/n 

5.84 

2.4% MORE 

1 

Radius, cm 
Rotor Exit Annulus 

22.86 

SAME 

• 

Area, cm^ 

656.7 

SAME 


Rotor Clearance, cm 

0.0 

SAME 

> 

Stator Solidity 

0.955 

SAME 


Rotor Solidity 

1.517 

SAME 


Stator Axial Chord, cm 

2.29 

SAME 


Rotor Axial Chord, cm 

2.79 

SAME 

1 

Stator TE Thickness, cm 

0.050 

SAME 


Rotor TE Thickness, cm 

0.060 

SAME 

Rotor Exit Swirl, deg 

-52 

SAME 


case represents only the first stage of a multi-stage turbine, so that 
this result need not be surprising. 

* Figures 17 and 18 portray the fixed geometry parametric study 
for rotor solidity and exit swirl, respectively. Unlike the previous 
case, predicted efficiency is strongly dependent on these parameters.. 

* In this instance, the original design efficiency and the measured test-rig 
efficiency were identical (86.0) as compared to 85.0 for the model prediction. 
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Figure 16: Case 4 Preliminary Geometry Optimization Study 
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Case 4 Efficiency vs. Rotor Solidity; Comparison of Code Pre 
diction. Original Design Value, and Test-Rig Measurement 
(Fixed Geometry) 






ROTOR EXIT SWIRL, 

DEGREES 

Figure 18: Case 4 Efficiency vs. Rotor Exit Swirl for Fixed Stator and 

Rotor Solidities (Fixed Geometry) 


Case 5 1 6 

Case 5 represents a compressor drive turbine originally intended 
for use in a compact automobile. The design is characterized by a large 
(49°) stator inlet flow angle, required to match the swirl distribution 
in the tangential entry inlet manifold. Pertinent model input parameters 
are given in Table 6. It should be noted that the model input corre- 
sponds to the so-called "smoothed and thinned profile" test (i.e., 
smoothed and thinned blades) of Reference 16. 
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Table 6; Case 5 Performance Model Input 


Parameter 

Design 

Equivalent 

Condition 

Test-Rig 

Condition 

Rotative Speed, RPM 

27673 

100% of design 

Inlet Temperature, °K 
Inlet Pressure, N/cm 2 

288.2 

320 

10.13 

8.0 

Specific Work, J/g 

44. A 

~5% less 

Mass Flow, kg/s 

0.325 

5% less 

Rotor Exit Mean 
Radius, cm 

5.00 

same 

Rotor Exit Annulus 
Area, cm 2 

36.47 

same 

Rotor Clearance, cm 

0.025 

same 

Stator Solidity 

0.557 

SAME 

Rotor Solidity 

1.803 

SAME 

Stator Axial Chord, cm 

1.17 

SAME 

Rotor Axial Chord, cm 

0.91 

SAME 

Stator TE Thickness, cm 

0.038 

SAME 

Rotor TE Thickness, cm 

0.038 

SAME 

Rotor Exit Swirl, deg 

-21.1 

SAME 


Figure 19 shows the geometry optimization portion of the para- 
metric study for Case 5. Predicted and actual turbine size are in good, 
but not excellent agreement. Figures 20 and 21 present the remainder of 
the parametric study. As with the previous case, efficiency appears to 
be sensitive to rotor exit swirl. Predicted efficiency at the test-rig 
conditions is 83. 5, which compares to a measured value of 82.5 and a 
design value of 85.0. 
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Figure 19: Case 5 Preliminary Geometry Optimization Study 



ROTOR SOLIDITY, 0^ 

Figure 20: Case 5 Efficiency vs. Rotor Solidity; Comparison of Code 

Prediction, Original Design Value, and Test-Rig Measurement 
(Fixed Geometry) 
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ROTOR EXIT SWIRL, 
DEGREES 

Figure 21: Case 5 Efficiency vs. Rotor Exit Swirl for Fixed Stator and 
Rotor Solidities {Fixed Geometry) 


Case 6 1 7 

Case 6 represents a single-stage "core 1 ' turbine for a turbofan 
engine and is characterized by relatively high hub-to-tip radius ratio 
and low aspect ratio. Due to relatively high requirements for work extrac- 
tion, Mach number levels are also high. The vane exit flow angle is flat, 
about 73° from axial; the vanes are untwisted and have a constant section 
profile. Pertinent model input parameters are given in Table 7. 

Figure 22 illustrates the geometry optimization calculation for 
Case 6. Predicted and actual turbine size are in excellent agreement. 
Figure 23 presents the rotor exit swirl parametric study for this case. 


Table 7: Case 6 Performance Model Input 


Design 

Parameter Equivalent Test-Rig 

Condition Condition 


Rotative Speed, RPM 
Inlet Temperature, °K 
Inlet Pressure, N/cm 2 
Specific Work, J/g 
Mass Flow, kg/s 
Rotor Exit Mean 
Radius, cm 
Rotor Exit Annulus 
Area, cm 2 

Rotor Clearance, cm 
Stator Solidity 
Rotor Solidity 
Stator Axial Chord, cm 
Rotor Axial Chord, cm 
Stator TE Thickness, Ci. 
Rotor TE Thickness, cm 
Rotor Exit Swirl, deg 


8081 

100% OF design 

288.2 

378 

10.13 

24.13 

76.84 

0.5% LESS 

3.708 

4% MORE 

23.5 

SAME 

562.5 

SAME 

0.030 

SAME 

0.929 

SAME 

1.487 

SAME 

3.81 

SAME 

3.43 

SAME 

0.127 

SAME 

0.127 

SAME 

-23.7 

-22.6 


Note that for the test-rig geometry the calculation model reaches the 
stator choke condition at a rotor exit swirl of approximately -27°. 

The measured swirl, on the other hand, is -22.6°, indicating that the 
flow is supersonic at the stator exit, a condition confirmed by the 
original design velocity diagrams. Strictly speaking, the model cannot 
be applied to this situation since it is constrained to stator exit Mach 
numbers less than unity. If one extrapolates the code prediction curve, 
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PREDICTED 



Figure 22: Case 6 Preliminary Geometry Optimization Study 


however, the model predicts an efficiency level very near the 88.6 
value actually measured. No solidity study was run for this case. 


Case 7 represents the first stage of a two-stage turbine designed 
to drive the compressor and' fan of a low-cost turbofan engine suitable 
for light aircraft application. The aerodynamic design is rather 
conservative, with low Mach number levels and low stator turning. The 
turbine was tested as both a two-stage and a single-stage unit. Data 
for the single-stage test are given in Table 8. 
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Figure 23: Case 6 Efficiency vs. Rotor Exit Swirl for Fixed Stator and 

Rotor Solidities (Fixed Geometry) 


Figure 24 indicates excellent agreement between predicted and 
actual turbine size, even though Case 7 is, like Case 4, the first 
stage of a multi-stage turbine. For fixed geometry input, Figure 25 
indicates a predicted efficiency of 91.1, compared to a test efficiency 
of 93.0. The original design value for this turbine was 87.0. Figure 26 
indicates that, within the exit swirl range studied, the predicted 
efficiency curve is virtually flat. 


Table 8: Case 7 Performance Model Input 


Parameter 

Design 

Equivalent 

Condition 

Test-Rig 

Condition 

Rotative Speed , RPM 

15336 

100 % OF DESIGN 

Inlet Temperature, °K 

288.2 

300 

Inlet Pressure, N/cm^ 

10.13 

13.79 

Specific Work, J/g 

45,83 

6 , 5 % more 

Mass Flow, kg/s 

1.989 

0 . 8 % MORE 

Rotor Exit Mean 



Radius, cm 

10.16 

SAME 

Rotor Exit Annulus 



Area, cm* 

233.5 

SAME 

Rotor Clearance, cm 

0.030 

SAME 

Stator Solidity 

1.049 

SAME 

Rotor Solidity 

1.469 

SAME 

Stator Axial Chord, cm 

1.91 

SAME 

Rotor Axial Chord, cm 

2.23 

SAME 

Stator TE Thickness, cm 

0.050 

SAME 

Rotor TE Thickness, cm 

0.050 

SAME 

Rotor Exit Swirl, deg 

- 26.1 

- 26.5 
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ROTOR EXIT 
MEAN RADIUS, cm 


24: Case 7 Preliminary Geometry Optimization Study 


TEST RIO 
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CODE 
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ORIGINAL 
DESIGN VALVE 


1.0 1.2 1.4 16 18 2.0 2.2 2.4 


ROTOR SOLIDITY, Or 

Figure 25: Casa 7 Efficiency vs. Rotor Solidity; Comparison of Code Pre- 
diction, Original Design Value, and Test-Rig Measurement 
(Fixed Geometry) 
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ROTOR EXIT SWIRL, 

DEGREES 

Figure 26: Case 7 Efficiency vs. Rotor Exit Swirl for Fixed Stator and 

Rotor Solidities (Fixed Geometry) 

Summary of Results 

Results for each of the seven cases considered in this study are 
summarized in Figure 27 and Table 9. 

Despite obvious limitations of the model such as the pitchline 
nature of the analysis, simplicity of the boundary layer assumptions, 
etc., the results for the seven cases detailed above indicate that the 
model is capable of providing preliminary aerodynamic performance data 
with an acceptable degree of confidence. Experience with the model 
has shown that it is important to set up the correct stator inlet flow 
conditions (Mach number and pre-swirl) in order to achieve optimum 
results. 
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MEASURED TOTAL - 
TOTAL EFFICIENCY 

Figure 27: Predicted vs. Measured Efficiency; All Gases 

Finally, a description of the model and code verification 
results have been presented at the 18th and 19th AIAA/SAE/ASME Joint 
Propulsion Conferences held in Cleveland, Ohio (1982) and Seattle, 
Washington (1983), respectively^ 9 ’ 20 . 

Appendices A-E contain supporting equations and analyses for the 
model. Appendix F describes the required model input; Appendix G con- 
tains a sample output; and Appendix H is a listing of the FORTRAN 


program. 


Table 9: Summary of Predicted vs. Measured Total- 
Total Efficiency; All Test Cases 


Case 

Predicted a 

Efficiency 

Measured 

Efficiency 

1 

91.4 

91.0 

2 

86.6 

87-88® 

3 

84.1 

83.2 

A 

85.0 

86.0 

5 

83.5 

. 82.5 

6 

89. 3 C 

88.6 

7 

91.2 

93.0 


a At measured rotor solidity and exit swirl. 
Depending on twisted vs. untwisted rotor 
configuration, 

Extrapolated from Figure 23 data, 
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APPENDIX A 

STANITZ METHOD SUPPORTING EQUATIONS 

The Stanltz method explicitly requires the following thermo- 
dynamic flow properties: 


t m 


a 

t; 


T" 

JM 

T'i * 

i 


M 


P" 

P M 


o i » D" 

p i p i 


Energy Equation: 


c p T M + 2gJ 


M_ _ 


Vi 

c p T i + 2gJ + W shaft 


( 1 ) 


Moment of Momentum: 


“shaft * gJ * (“W 


( 2 ) 


Upon substitution of Eq. (2) into Eq. (1), we hav< 


t m_. < u Vi - < u Vm 

t ; 9 jt P T i 


(3) 


which can be written, after some manipulation, as 


T M 


= 1-2 


XjlI 


Y + 1 


l a crj 


ajL, 
crj 


- r m 


rv 


a 

L crj 



M 


T i 


(4) 


Eq. (4) can be combined with Eq. (5) of the main text to yield 


I 
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(5) 


Eqs. (1), (2), and the law of cosines relationship between absolute 
velocity, V, and relative velocity, W, can be combined to yield 






(1 - RJ 2 ) 


( 6 ) 


From the definitions of absolute and relative total temperature we can 
write 


with 


T" 

M 


U M 

T* * - M 

M 2gJc r 




(7) 


( 8 ) 


Eqs. (7) and (8) may be combined to give 
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Substituting Eq. (10) Into Eq. (13), we have 
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hr 




act 




act 
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(14) 


Equation (14) Is usable In both rotating and non-rotating reference 
frames. For stators, Eq. (14) reduces to the familiar assumption of 
linear total pressure loss, I.e., 



+ 



(15) 


The quantity (Pg/Pp act I s determined from the Stewart mixing hypothesis 
(Appendix B). 

Finally, since 



P M 

*T 


p ? 



we have 



( 16 ) 


where 
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III 

Y+l 


r wj 

a »l 

a cr 


2 

M 


must be evaluated for both the suction surface and the pressure surface 
of each blade. 


APPENDIX B 
PROFILE LOSS 
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Profile loss is defined here as a combination of frictional 
effects arising from the flow of a viscous fluid over a solid surface 
and the subsequent downstream mixing of the suction-surface and the 
pressure-surface boundary layers. Pressure drag, which results from 
the flow of fluid past a finite-thickness trailing edge, is implicitly 
contained in the Stewart analyses (Appendix C). 


Laminar Regime 

Laminar boundary Uyer properties are calculated as follows: 

1. The freestream velocity function V(x) and its de'ivative dV/dx are 
known . 

2. The momentum thickness, 9 (x) is calculated by 


0(x) _ X 
A ‘ V„ 


-3 



1/2 


d(f) 


1 


0 ) 


where - free stream velocity at the cascade exit, i.e., at x = s, 
C f = 1.328 (Re^)"^ 2 ; Re^ = Reynolds No. 
x t = location where boundary layer transition occurs 

3. The parameter Z, given by 


Z = 



v 


(2) 


where v s dynamic viscosity, is defined. 
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4. The parameter K, given by 



is defined. 

5. The shape factor A, given by 

k - (JL J_ a 1 a2\2 . 

K ‘ '315 “ 945 A “ 9072 A ^ A 
is obtained by iteration. 

6. The displacement thickness, &*, is then calculated from 

_L _ _!_ A 
6* _ 10 120 A 
e " JT 1 7 ~ 1 ,2 

315 * 945 A “ 9072 A 


(3) 


(4) 


(5) 


Transition 

Reference 6 indicates that the point of instability for boundary 
layers in a pressure gradient (that is, the point at which disturbances 
begin to amplify) can be determined through consideration of a "critical 1 ' 
Reynolds number, Re cr , based on boundary layer displacement thickness. 
Further, this critical Reynolds number is a function of the shape factor 
A. For the present model, two assumptions are made: 

(1) Re cr (A) is obtained from a curve fit of data given in [6]. 

(2) The point of instability is assumed to coincide with the 
point of transition. 
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Turbulent Regime 

Boundary layer momentum thickness is computed from an equation 
appearing in References (3) and (4): 

e(x) = 


0.231 


f L 


a* 


V 

< 

cr' 


cr 


(1+HT 


1.268 


ipvj 


0.268 


(l_A)°* 467 dx| 


^ cr' 1 cr' 


(1+H b ) 


10 


0.6T8(2n+l)' 


where 0 = boundary layer momentum thickness at the blade trailing edge 
u = viscosity 

n = exponent in the boundary layer power-law velocity profile, 
taken as 1/7 

x = distance along blade surface 
The quantity "A" is defined as 


A 


111 

Y+l 


cr 


(7) 


and "H" is the boundary layer form factor, defined as 


V (2m+l )A m 
m =0 \2m+1 )n+l 

oo 

l Q [{2m+l )n+l ][2(m+l )n+lj 


( 8 ) 


.7886 

( 6 ) 


from Reference (7). All other parameters used in Equation (6) are 
defined elsewhere. 
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ft 

Since the form factor is also defined as 

» 

H * T < 9 > 

I where. 6 * is boundary layer displacement thickness. Equations (6), (7), 

{8} and (9) adequately describe all the pertinent turbulent boundary 
layer parameters required to calculate profile loss. Profile loss 
s itself is calculated from the Stewart Mixing Hypothesis of Reference 

(4), which is detailed in Appendix C. 

ft 

ft 

I 


i 




I 


I 



APPENDIX C 

STEWART MIXING HYPOTHESIS* 4 * 

Following Stewart, we define stations at the cascade inlet, 

1 ) cascade exit, and (somewhere) downstream of the cascade exit, as shown 

in Figure C-l. 


) 



Figure C-l 


Continuity: 

1 

cos a-, j(pV) 1 d(~) = ( pV ) 2 cos a 2 (1) 

0 

Axial Momentum: 

1 

gP 1 + cos 2 oJ (pV 2 ) 1 d(^') = gP 2 + cos 2 a 2 ( P V 2 ) 2 (2) 

0 

Tangential Momentum: 

1 

f ? ii p 

sin a 1 cos aJ (pV )^d(—) = sin a 2 cos a 2 (pV ) 2 (3) 

0 
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Define 



(4-a) 


(4-b) 


where the subscript "fs" refers to freestream conditions (assumed 
isentropic). 

We shall define a set of reference (total) conditions at station 1 
and assume adiabatic flow throughout. Eqs. (1) and (4-a) can then be 
written as 


(1 - 6 * - 


V 


f pV ] 

♦ 

[ pv f] 

P 2 

Ki 

fs»l 

Krl 

2 P fs,l 


where the subscript "x" refers to the axial direction. 
Eq. (4-b) can be combined with (4-a) to become 


1 


9 * = (1 - <$* - 6 te ) - 


f i 

f > 

_e__ 


r v ] 

o 

Pfs^ 

i 



d(f) 


(5) 


( 6 ) 


Eq. (2) can then be written, after algebraic manipulation, as 





while Eq. (3) can be written as 
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(!-<$*- 6 te - 0*) 


_ V V ] 


„ V V.,1 


r ?' 

P x u 


P x u 


2 

p 1 al~ ITT 

cr cr 

fs,l 

p a cr a cr 

2 

p* 
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(8) 


where the subscript "u" refers to the tangential direction. We shall 
define 

,2 


Y-l 

r v 

Y+l 

a cr 

V 


fs,'l 


(9) 


and set 


J 1 ’1 


P fs,l p fs,l T fs,l 


v] (1 ' A fsl } 
t p J fSfl TS,I 


Now solve 
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Eq. (5) for p 
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00) 


P 2 

2 .= (1 - 5 * - S ) 
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01) 


Substitute Eqs. (11) and (10) into Eq. (7) and solve for 


fv x l 
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( 12 ) 


where 
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(1 * A fs,l> (l?) * 0 * S * ' S te • 9 *> 


C = 


f V > 

a 1 " 

l cr ) 


«te> 


cr 


fill 


fs.l 


03) 


Substitute Eq. (11) into Eq. (8) and solve for 


fv. 
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cr 

V / 


r V 
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cr 


(1 - 6 * - 6 te - 0 *) 
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(’-«*- 4 te ) 


15 . 1=0 


(14) 


Eq. (12) can then be written as 
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which has the solution 
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For stators , Pjr sl = Pg so that 
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fs.l 


1 a i 
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(19) 


For rotors , we will assume no radius change between the cascade exit 
(station 1) and the mixed plane (station 2); there may, however, be a 
change in radius between the cascade inlet (station 0) and the cascade 
exit. We have 
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where -gn is an isentropic change given by 
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Thus, we have 
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mix 


where the subscript "mix" refers to station 2 conditions. Note that 
the LHS of Eq. (20) appears In the efficiency equation, Eq. (11) of 
Appendix E. Pressure loss terms are thus functions of <5*. and 9*; no 
Isentroplc terms appear. 




APPENDIX D 

BLOCKAGE CALCULATIONS 
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The effect of leading edge and trailing edge blockage Is deter- 
mined from the continuity equation and the assumption that moving from 
blocked to unblocked flow (and vice versa) leaves the tangential momentum 
(velocity) unchanged. Thus, one can easily show that, for the trailing 
edge. 
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and for the leading edge 
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where 0j^, 0^, 3 M j x = trailing edge blade angle, leading edge blade 


angle, and appropriate dowirstreom ouiu upstream 


flow angles, respectively 
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tjE» = trailing edge and leading edge blade thickness, 
respectively 
s - blade spacing 

w 

The critical velocity ratio, tt— , is always defined relative to the 

a cr 

blade, whether for stators or for rotors. 

It should be noted that the present model uses Equations (1) 
and (2) to modify (assumed) blade loadings to account for finite leading 
and trailing edge blade thickness. 
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APPENDIX E 


STAGE EFFICIENCY 



By definition, the adiabatic efficiency of the turbine 
expansion process shown schematically in Figure E-l is 


T 1 - T 3A 


n T- T ' TJ - T£ 


which can be written as 
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Denoting entropy changes by As, we have 


As total As stator + As rotor 
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As rotor = as 2%3*^3J = V" 
From Eq. (7) we can write 
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so that 
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Thus, combining Eqs. (3), (5), (6), and (9) we have 
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Substituting Eq. (10) into Eq. (2) yields the stage efficiency: 
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APPENDIX F 
MODEL INPUT 


The calculation model is set up to use the English system of 
> units (ft, sec, lb m ). Required program input is as follows: 



1 . 

GAMMA: 

ratio of gas specific heats, assumed constant throughout 
the calculation 

> 

2. 

RGAS: 

gas constant, ft*lbs/lb m /°R 


3. 

SPEED: 

rotor rotational speed, RPM 


4 . 

D 1P: 

turbine inlet total pressure, psia 

) 

5. 

TIP: 

turbine inlet total temperature, °R 


6. 

RWORK: 

turbine work output required, BTU/lb m 


7. 

WFLOW : 

mass flow rate, lb m /sec 

m 

> 

8. 

RM1 , RM2, RM3: pitch line radius at stator inlet, rotor inlet, 

and rotor exit, respectively (inches). If DSTRES = 0, 
input RM as zero also. 

) 

9. 

AREA! , 

AREA2, AREA3: annular flowpath area at stator inlet, rotor 
inlet, and rotor exit, respectively (sq. inches). If 
DSTRES - 0, input AREA as zero also. 

> 

10. 

CLEAR: 

rotor tip clearance, inches 


11. 

TETS: 

stator trailing edge thickness, inches 


12. 

TETR: 

rotor trailing edge thickness, inches 

) 

13. 

ALP2: 

stator exit flow angle, degrees. IF ALP2 = 0, angle will 
be calculated assuming zero rotor exit swirl. See note on 
ALP3. 

) 

14. 

ZWFS: 

stator Zweifel coefficient. Assumed 0.8 if none is input. 


15. 

ZWFR: 

rotor Zweifel coefficient. Assumed 0.8 if none is input. 


> 
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16. 


17. 


18. 


19. 


20 . 

21 . 


22 . 


23. 


24. 


ALP1 : stator inlet flow angle, degrees 

CH0R0S, CHORDR: stator and rotor axial chord, respectively (inches). 

If input as zero, chords will be calculated by aerodynamic 
optimization. 

TEIS, TELR: stator and rotor leading edge thickness, respectively 

(inches) 

OSTRES: maximum allowable disk stress, psi. If DSTRES f 0, program 
will calculate a range of possible turbine sizes and 
performances. If DSTRES = 0, program will perform an 
aerodynamic optimization for a single case only; RM and 
AREA must be input in this case. 

VEXIT: rotor exit axial critical velocity ratio, V /a' 

X CP 

ALP3: rotor exit swirl angle, degrees. ALP3 and ALP2 cannot 

both be input; if both are, only ALP3 is used. If neither 
ALP3 or ALP2 is input, program will calculate ALP2 assuming 
ALP3 = 0. 

IPLOT: key for blade plot generation. Set IPLOT = 0 for no plot. 

Note: user will have to tailor program for local plot 

subroutines . 

RHODI: disk material density, lb m /ft 3 

RHOBL: blade material density, lb m /ft 3 

blade-disk attachment region material density, lb m /ft 3 
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PROGRAM TO COMPUTE ADIABATIC EFFICIENCY FOR 
SINGLE STAGE AXIAe GAS TURBINES 

WORK PERFORMED UNDER NASA GRANT NO. NSG 3293 
NASA LEWIS RESEARCH CENTER 
CLEVELAND* OHIO 

DIMENSION URELS<S1 ) x VR£LP(5l ) ,AM(51 ) ,RMSTR(51 > ,HMSTR(31 ) ,TMT2P(5l ) 
DIMENSION PMPIP <51 ) xGRAN< 51 ) xSIG(Z) ,CORD(2) xTMT2(5l ) xVUCRM(51) 
DIMENSION WUCRM(31 ) xWCRM(51 ) .BETAMOl ) x£PSM(51 > xTPTPM(5l ) xTHIK(5l ) 
DIMENSION THIKR(31 > x THIKL<51 > rTHETA<31 ) x RTHET (51 ) xR(5l ) xDIFFl (2) 
DIMENSION HMS(51) xHMP(51 ) ,WS(51 ) xWP<51 > , DTE ( 2 ) x PLOSS ( 2 ) xDIFF2(Z) 
DIMENSION YPRO<2) »YSF(2) x YT0T(2) xCHRD(200) , YT0TL(200) xREYNQ(200> 
DIMENSION FSL (51 ) x FPL (51 ) ,BNUS(51 > xSNUP(51 ) xDUDXS(51 ) ,DUDXP(5l ) 
DIMENSION FS(31) rFP(51 ) »GS (51) rGP( 31 ) xSEFFYUl ) xR0S(51 ) ,ROP(5l ) 
DIMENSION CXLEN(51 ) ,RADEX(2l > ,EADIAB(2l ) .FAD I AS (200) .DPS I (2l ) 
DIMENSION RADEXX(21 ) x TITLE (20) x VDP(200> 

4000 READ ( 45 , * x END -999 ) GAMMA r RGAS x SPEED x PlP » T1 P x RWORK » WFLOW 
READ(43.*)RM1 » RM2 » RM3 » AREA 1 » AREA2 x AREA3 x CLEAR 
READ ( 45 x # ) TETS . TETR , ALP2 x ZWFS x ZWFR x ALP 1 x CHORDS 
READ ( 45 x * ) CHORDR x TELS » TELRx DSTRES x VEX I T , ALP3 x I PLOT 
READ ( 45 » * ) RHOD I x RHOBL , RHO AT 
WRITE(6r210l ) 

WRITE(6x2102) 

WR I TE ( 6 x 2 1 03 ) GAMMA . RGAS , S PEED x WFLOW , P 1 P , T 1 P . RWORK 
WRITE(6r2104> 

WR I TE ( 6 x 2 1 05 ) RM 1 x RM2 x RM3 x ARE A 1 x AREA2 x ARE A3 , TETS x TETR x CLE AR 
WRITE ( 6 x 2201 ) 

WR ITE ( 6 r 2202 ) ALPl x ALP2 x ZWFS , ZWFR x CHORDS . CHORDR 
WR I TE ( 6 x 2204 ) TELS x TELR x DSTRES , VEX I T , I PLOT x Ai.P3 
WR I TE ( 6 » 2245 ) RHOD I x RHOBL , RHOAT 
WRITE (6x2203) 

IF( DSTRES. Efl. 0.0) WRITE(6x2242) 

IF (DSTRES .EG. 0.0) WRITE ( 6 x 2243 ) 

GAMM1»( GAMMA- 1 . >/(GAMMA+l . ) 

GAMM2*( GAMMA-1 .) /GAMMA 
GAMM3 ■ 1 . / ( GAMMA- 1 . ) 

SETAL3»ALP3*3. 14159/ISO. 

IF ( DSTRES • NE .0.0) RMl-0.0 

I F ( DSTRES . NE . 0 . 0 ) RM2=0.0 

IF (DSTRES. ME. 0.0) RM3=0.0 

IF ( DSTRES. NE. 0,0) AREA1=0.0 

I F ( DSTRES . NE . 0 . 0 ) AR£A2=0.0 

IF ( DSTRES. NE. 0.0) ARE A3* 0.0 

IF (DSTRES. NE. 0.0) ALP1*AlP1*3. 14159/180 . 

ALP3=At.P3*3. 14159/180. 

KPAS5*0 

RBEST=0.0 

190 DO 19 1*1x200 
CHRD( I ) *0. 

YTOTL( I ) *0 . 

REYNO ( I ) *0 . 

19 CONTINUE 
1PRINT*0 
NPASS*1 

IF (DSTRES. EG .0.0) GO TO 175 

1 80 CALL FD I SK ( GAMMA x VEX I T , GAMM 1 x GAMM2 x GAMM3 x NPASS x K PASS x RGAS x SPEED x 
1 P 1 P , T1 P , RWORK x WFLOW x ALPl x RM1 x RM2 x RM3 , HI r H2 x H3 x T2P x U2 » UACR2 » U3 x 
2V 1 x VU 1 , ETA x P2P x T3P x P3P x PRTURB ? At.P2 x V2 x VU2 x VWRK2 , VWRK3 x UACR3 x VU3 x 
3TPTP2 x T2PP x PDPP2 x T32PP x T3PP x WU2 x WUPP2 x 8ETA2 x WPP2 x WU3 x WUPP3 x V3 x 
4ALP3 r VX3 x WXPP3 x WPP3 x BET A3 , P32PP x PLOS 1 , PLGS2 x RH3 x RT3 , RTS x ACR 1 x ADP3 x 
5H 1 A x H2A x H3A x VX2 x WFACTM x WFACTH x UT3 x STG ACC , RBEST x T2PPH x ACR3 ) 

... GO TO 11 
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1 
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1 
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10 
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» 


> 


210 


> 


> 


i 


* 212 
213 


BLP2-ALP2 

RMl-RHl/12. 

RM2-RM2/12. 

RM3-RM3/12. 

Hl»<AREAl/l44.)/(2.*3.l4l59*RMl ) 

H2- < AREA2/ 144 . ) / < 2 . *3 . 1 4l59*RM2 ) 

H3-( AREAS/ 144. )/ (2 .*3. 14159*RM3> 

H1A-H1 

H2A-H2 

H3A-H3 

T2P-T1P 

TSTD-518.7 ' 

ASTD* ACR I T ( GAMMA * RGAS • TSTD ) 

RHSTD*2ll6.22/(RGAS#TSTD> 

CP- < 1 . /GAflM2 ) *RGAS/778 . 26 
U2*3. l4159*SPEED*RM2/30. 

ACR 1 “ACR I T ( GAMMA » RGAS » Tl P ) 

UACR2=U2/ACRl 
U3»U2*RM3/RM2 
UT3-U3* ( RM3+H3/2 . ) /RM3 
ALP 1 -ALP 1*3. 14159/180. 

BLPl=ASS<ALPl ) 

01 -VELJT ( WFLOW » Tl P r GAMM1 r GAMM3 r BLPl 7 AREAl 7 P 1 P 7 RHSTD r ASTD ) 

VUl a Vl *S I N ( Ai_Pi ) 

SWPM* < U2**2 ) / ( 25036 . 62*RNGRK > 

W0RKC=1./SWPM 

ETA a (0. 92*SWPM> / ( E 0.0227) 

P2P-0. 93*P1 P 
T'spaTi p-ounpk /rp 

P3P-P1P*< <1.-<1.-T3P/T1P)/ETA)*#<1./GAMM2) ) 

PRTURB-P1P/P3P 

KALP“0 

IF ( 8LP2 . EG .0.0) GO TO 210 
ALP2-8LP2 

BNGCHK a WFLGW*SGRT ( Tl P/TSTD > / ( ( < 1 . -GAUM 1 ) **GAMM3 ) *AREA2/ 1 44 . *P2P/ 
114. 696*RHSTD*ASTD ) 

ANGCHK 3 AC05 < BNGCHK ) 

ANGCHK=ANGCHK*l80./3.l4lS9 

IF( IPRINT.EG. 1 ) WRITE ( 6 . 2210 ) ANGCHK 

IF1ALP2.GT. ANGCHK. AND. IPRINT. EG. 1) WR I TE < 6 r 2200 ) ANGCHK 
Al_P2*ALP2*3. 14159/180. 

ANGCHK 3 ANGCHK*3. 14159/.180. 

IF (ALP2.GT. ANGCHK) ALP2=ANGCHK 

V2 3 VELIT ( WFLOW , Tl P , GAMMl , GAMM3 » ALP2 7 AREA2 r P2P , RHSTD , ASTD ) 
VU2 a V2*S IN ( ALPZ ) 

VWRK2»VU2*ACR1 

VWRK3»RM2/RM3*VWRK2-32 , 174*778 . 26*CP*T1 P/V2* ( 1 . -T3P/T1 P ) 

I F ( V WRK 3 . NE .0.0) KAl_F*l 
GO TO 220 

I F ( SET AL3 . EG . 0 . 0 ) GO TO 211 
AL2 3 A6S ( SETAL3 ) 

V3“ VEL I T < WFLOW »T3P. GAMM 1 » GAMM3 7 AL3 r AREA3 1 P3P» RHSTD , ASTD ) 

VU3 » V3*S I N ( SETAL3 ) 

VX3»V3*C0S ( SETAL3 ) 

ALP3-SETAL3 

ACR3» ACR I T ( GAMMA , RGAS t T3P ) 

VWRK3*VU3*ACR3 

VWRK2-32 . 1 74*778 . 26*CP*T1 P/U2* ( 1 . -T3P/T 1 P ) +RM3/RM2*VWRK3 

VU2-VWRK2/ACR1 

IF( VU2. GT .1.0) GO TO 212 

ALP2-ANG I T ( WFLOW » Tl P r GAMM 1 1 GAMM3 > VU2 7 AREA2 7 P2P 7 RHSTD r ASTD 7 0 ) 

CHEK 3 VU2/ S I N < ALP2 ) 

IF<CH£K.LT .1,0) GO TO 213 
WRITE(672239) 

GO TO 4000 
V2-CHEK 


VX2“V2*C0S < AuP2 ) 

UACR3-U3/ACR3 

GO TO 214 87 

UNRK3-0.0 

VNRK2*32. 174*778. 26*CP*T1 P/U2* < 1 . -T3P/T1 P ) +RM3/RM2*VkRK3 

UU2»UWRK2/ACRl 

IF( VU2.GE. 1.0) GO TO 15 

AlP2»AnGIT < INFLOW r 71 P , GAMMl » GAMH3 , UU2 , AREA2 .P2P, RHSTD , ASTD , 0 ) 
CHEK«VU2/«3IN<ALP2) 

IF(CHEK.LT.l.O) GO TO 12 
VWRK3-UWRK3-10. 

KAuP“l 
GO TO 13 
V2-C.HEK 

ACR3 « ACR I T ( GAMMA r RGAS » 73 P ) 

VX2»U2*C0S<ALP2> 

UACR3-U3/ACR3 

VU3«VWRK3/ACR3 

TPTP2"1 . -GAMMl #(UACR 2**2. )*<2 .*VU2/UACR2-1 . ) 

T2PP»7PTP2*T2P 
PDPP2»TPTP2**< 1./GAMM2) 

T32PP-1 . -GAMMl *<UACR2**2. )*<1 ,/TPTP2)#< i .-( i RM3/RM2 > **2 . ) ) 

T3PP*T32PP*T2PP 

ADP2-ACR IT < GAMMA , RGAS , T2PP ) 

ADP3»ADP2*SGRT(T32PP> 
kU2 “VWRK2-U2 
WUPP2-WU2/ADP2 

B£7A2»ASIN< ( VU2-UACR2 > / ( SGRT ( V2**2+bACR2**2-2 . * VU2*UAC-R2 ) ) ) 
WPP2»WUPP2/SIN(BETA2) 

WU3*VWRK3-U3 
WUPP3*WU3/ADP3 
I F ( SET AL3 . NE .0.0) GO TO 216 
IF< K ALP. EG . 1 ) GO TO 14 

V3*UELIT ( WFLGN , T3P » GAMM 1 , GAMM3 r 0 . 0 . AREAS , P2P , RHSTD . ASTD ) 

ALP3=0.0 
GO TO 17 
UUM=UU3 
VUM»A8S(VUM) 

ALP3 3 ANGIT ( WFLQW» T3P , GAMMl , GAMM3 , VUM , AREAS . PSP, RhS . ASTD , 1 > 

IF< VU3. LT .0.0) ALP3=-ALP-3 
V3 a VU3/SIN< ALPS ) 

UX3 a V3*C0S < ALPS ) 
wXPP3*VX3*ACR3/ADP3 
WPP3=SGRT ( WXPP3**2+wUPP3**2 ) 

WFACTM* 25036 . 62*RN0RK / ( U3**2 ) 

.WFACTH-WFACTM* ( ( RM3/ ( RM3-H3/2 . ) ) **2 > 

BETA3a AS IN { WUPP3/WPP3 ) 

IF ( NPASS . NE. 1 ) GO TO 11 
P32PP a 732PP**< 1 . /GAMM2 ) 

PL0S1 =0.98 

PLQS2 a P32PP 

CONTINUE 

IF ( I PRINT .EG. 1 ) WRITE(6,2100) 

IF(RBEST.GT.O.O. AND. I PRINT. EG. 1 > NRITE(6r2233) 

IF< I PRINT . EG . 1 ) HRITE (6*2106) 

IF< I PRINT . EG . 1 ) WR I TE ( 6 * 2 1 07 ) VI * 0X3 » VU 1 » VUG * V2 * V3 * VX2 , HPP2 

IF < I PRINT . EG . 1 > HRITE < 6 * 21 08 ) VU2 * k*PP 3 , P2P, UACR2 , PLOS 1 * UACR3 * PL0S2 * 

1 P2P. T3P , UT3 

IF(DSTRES.£G.O.O.AND.IPRINT.EG. 1 )HRITE<6*2097> 

IF ( DSTRES . EG . 0 , 0 . AND . I PRINT . EG . 1 > HRI7E < 6 * 2240 ) NFACTM , wFACTH , PRTURB 
IF (DSTRES.EG.O.O) GO TO 2lS 
I F < RBEST . EG . 0 . 0 ) RADEX ( K PASS+1 ) =RM3 
RH3 a RH3*l2. 

RH2=RH3 

RH1*RH3 

RT3»RT3*l2. 
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ANGZ2*SETA3 

I AnGGAM* ABS ( BETAZ+8ETA3 ) /2 . 

ZI*£IFL*0.30 

IF(ZWFR.NE.O.O) ZWEIFL=ZkFR 
HEIGHT* <H2A+H3A) *6. 

CHORDB=CHORDR 
R3STR*RM3/RM2 
I H3STR»H3/H2 

VUCRI *VU2 
VUCR0=VU3 
AIN*ALP2 
VCRI*V2 
RIN»RM2*12. 

TIPDP-1 ./TPTP2 
' P0PI-PLGS2 

| T0TIA*T3P/TI P 

T0TIR-T32PP 

) PIPDP*1 . /PDPP2 

TPTPM(1)*TPTP2 
/ WCRM<l)*i«PP2 

> TIVIS*T2PP 

.1 RH02*l44.*PDPP2*PLQSl*PlP/(RGAS*T2PP> 

CALL VI SCO ( T3PP > VISRE ) 

] RH0RP*RHO2*P32PP 

! RHORE*RHOV ( GAMMA , WPP3 ) /NPP3*RH0RP 

i VELRE=HPP3*ADP3 

) VELSL-HPP3 

UELSL*wUPP3 
TL£*TELR 
VELBL0=WPP2 
BE7ALE*B£TA2 
TTE°TETR 

| BETATE=BE7A3 

60 CONTINUE 

IF ( DSTRES. £Q .0. 0 )GO TO 251 

IF(CHORDB.EG.O.O.AnD.RBEST.GT.O.O) GO TO 250 
IF ( CHORDS . NE . 0 . 0 . AND . RBEST . GT . 0 . 0 ) GO TO 252 

S0L£ST=2 . *COS ( ANGZ2 ) *S I N ( ABS ( AnGZI -AnGZ2 ) ) / ( ZHEIFL*COS ( AnGZI ) ) 

> CH0RD=0 . 26S07*H£IGHT*S0LES7/CQS < AnGGAM > 

ITRIG* 1 

4 I CHORD* 1 

GO TO 500 

251 I F ( CHORDB . EG .0.0) GO TO 250 

252 CHGRD=CHORDB 

v ITRIG* 1 

I CHORD =1 
GO TO 500 
250 I CHORD* 1 

ITRIG=0 

CHORD* HEIGHT/ 4. 

CHMAX=HEIGHT*3.0 
' 500 KZWFL»0 

119 S0LID*2.*C0S <ANGZZ)*SIN(A8S< ANGZI -ANGZ2) >/<Zk£lFL*COS (ANGZI ) ) 
SPAC=CHORD/SOLID 

CALL eLOCK ( GAMMA , TTE , SP AC , BETATE , UELBL , BETTE * VTE ) 

BETLE*BETALE 

VL£*VELBLO 

> EXTRM*0. 3 

CALL DRANGE < GAMMA > VLE , VTE , EXTRM , DST ART , DCHNGE ) 

DPBLD*DSTART 
ASLD*DCHNGE 
DELBD=DCHNGE 
K* 1 

- •* KK =0 

117 CONTINUE 

118 IF(IBLD.£G.2> GO TO 115 ... _ . 
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RT2=RTS*l2. 

RT1 ®RT2 

AREAR»3. 14159*<RT3**2-RH3**2> 

AREAS "3 . 14 159#< RT2**2-RH2**2 ) 

IF< I PRINT . EG. 1 ) HRITEl 6. 2109 ) RH3 .RT3 r RH2 . RT2 » RHl » RTl » AREAR > AREAS » 
1STGACC 

I F i I PR I NT . EG . 1 ) HR I TE < 6 , 2 1 46 ) HFACTM , kFACTh , FRTUR6 
215 ALP1 *AlP1*180./3. 14159 

ALP2=ALP2*l30./3. 14159 
BETA2=BETA2*180./3. 14159 
BETA3»BETA3*l80./3. 14159 
ALP3»AuP3*l80. /3. 14159 
IF( IPRINT.EG. 1 ) kRITE<6r2098> 

IF( IPRINT.EG. 1 ) kRITE<6,2l41 ) ALP1 . AlP 2 . BETA2 » BET A3 » At.P3 
ALP1 =AlP1*3. 14159/180. 

ALPZ»ALP2*3. 14159/180. 

SETA2»BETA2*3 .14159/180. 

BETA3=BE7A3*3 .141 59/ 1 SO . 

ALP3«ALP3*3. 14159/180. 

IF< IPRINT.EG. 1 ) kRITE<6,2100> 

COMPUTATION OF TRIAL VELOCITY DIAGRAMS IS COMPLETE 
!BLD*1 

20 IF( IBLD.EG.2) GO TO 50 

IF< IPRINT.EG. 1 ) kR I TE(6»2llO) 

RPM-0.0 
AnGZ1»AlP1 
ANGZ2*ALP2 

AnGGAM“ABS ( ALP l +ALP2 ) /2 . 

Zk£IFL=0.8 

IF(ZWFS.NE.O.O) ZHEIFL=ZkFS 
HEIGHT* (H1A+H2A) #6. 

CHORDS “CHORDS 
R3STR*RM2/RM1 
H3STR=H2/Hl 
VUCRI=VUl 
VUCR0-VU2 
AIN=Al_Pl 
VCRI =V1 
RIN*RM1*12. 

TIPDP° 1 .0 
P0PI=PL0S1 
TOTI A=1 .0 
TOTIR* 1 . 0 
PIPDP*1 .0 
TPTPM ( 1 >*1 .0 
kCRM(l)=Vl 
TIVIS»71P 

RH02 = 144 . #P1 P/ < RGAS#T1 f-*) 

CALL V I SCO < 7 1 P * VI SRE ) 

RH0RP*RH02*P2P/P1 P 
RHORE=RHOV ( GAMMA - 02 ) /V2*RH0RP 
VELRE * V2*ACR 1 
VELBLaMZ 
UEL8L*VUCR0 
TLE-TELS 
VELBL0=V1 
BE7ALE*ALP1 

IF< ALP1 . EG .0.0) BETAuE=0. 00001 
TTE*7ETS 
BETATE=ALPZ 
GO TO 60 
50 RPM*UACR2 

IF< IPRINT.EG. 1 > HRITE(6r2111) 

ANGZ1 *BETA2 


35E5! 
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CALL VNEG ( VLE p VTE p EXTRM r DPBLD r VRATIO) 

DO 21 1=1 p51 
AI*( 1-1 )/50. 

21 CALL D 1ST < VLE , VTE r DPBLD , AI r VRELS (Dp VRELP < I > p DUDXS (Dp DUDXP (I)p 
1 EXTRM p VRATIO) 

GO TO 116 

115 CAlL VNEG( VLE » VTE p EXTRM , DPBLD p VRATIO ) 

DO 22 j=1 p51 
Aj»( J-D/50. 

CALL DIST(VLEpVTEpDPBLDpAjpVRELS( J) ,VRELP( J) pDUDX5( J) pDUDXP( J) p 
IEXTRMpVRATIO) 

22 CONTINUE 
116 DO 30 J=1 p51 

AJ= J-l 

Ad( J ) “AJ/50. 

RMSTR ( J ) = < R3STR- 1.0) *AM ( J ) + 1 . 0 
HMSTR ( J ) “ ( H3STR- 1 . 0 ) *AM ( J ) + 1 . 0 

TMT2P ( J ) = 1 . O-GAMMl# (RPM**2 ) *TI PDP* ( 1 . - ( RMSTR ( J ) **2 ) ) 

PMP I P ( J ) *TMT2P( J ) ** ( 1 . /GAMM2 ) + ( PGP I -TOT I R** ( 1 . /GAMM2 ) ) *Art ( J ) 
PMP2S =PMPI P ( J ) # ( ( 1 . /TI PDP# ( 1 . O-GAMMl# ( VRELS ( J ) **Z ) ) ) 

1**( 1 . /GAMM2 ) ) 

PMP2P =PMPIP(J)*((1. /TI PDP* ( 1 . O-GAMMl * ( VRELP ( J ) **2 ) ) ) 

1**( 1 . /GAMM2) ) 

IF( IBLD . EG . 2 ) GO TO 31 
DELP=PMP2P-PMP2S 
GO TO 32 

31 DELP=PMP2S-PMP2P 

32 GRAN C J ) =HMSTR ( J ) *RMSTR( J ) *DELP 

30 CONTINUE 

SUM=S I MP2 ( GRAN , 0 . 0 , 0 . 02 p 50 ) 

Si I Gil =(2.*GAMMA/(GAMMA+1 .0) ) *RHGV < GAMMA p VCR D *C0S < AIN )*( R3STR* 

1 VUCRG*SGRT ( TOT I A ) -VUCR I) / SUM 
SIG(K)=SIGM 
160 DPBLD =DPBLD+DELBD 
K = K + 1 

IF ( K . EG. 2 ) GO TO 117 
DPBLD=DPBLD-2. *D£LBD 
K = 1 

u u 3 kk 4. i 

IF(KK . GT. 20 ) GO TO 125 
DIFF1 ( 1 )*ABS(SQLID-SIG< 1 ) ) 

DIFF1 (2)=ABS(S0LID-SIG(2) ) 

DIFF2(1)=S0LID-SIG(1) 

DIFF2( 2 ) =SOLID-SIG (2) 

IF(DIFFKl).LE.O.Ol) GO TO 120 
IF(DIFF1(2).LE.0.01> GO TO 121 
GBLD*DIFF2< 1 )*DJFF2(2) 

IF(GBLD.LT.O.O) GO TO 122 
IF(DIFF1 (2) .GT.DIFFKl ) ) GO TO 123 
DPBLD=DPBLD+ABLD*1 .0 
GO TO 117 

122 DELBD=DELBD/2. 

ABLD-ASLD/2. 

GO TO 117 

123 DELBD=-DELBD 
ABLD=-A6LD 
GO TO 117 

120 SI GMA=SIG ( 1 ) 

IF( IBLD. EG. 2) GO TO 126 

CALL VNEG ( VLE p VTE p EXTRM , DPBLD p VRATIO ) 

DO 127 1=1 p 51 
AI= ( 1-1 ) /50. 

127 CALL DIST( VLE pVTEp DPBLD r AI p VRELS ( I ) p VRELP ( I ) r DUDXS ( I ) p DUDXP ( I ) p 
IEXTRMpVRATIO) 

GO TO 305 

. 126. CALL .VNEG.(VLEpVTE,EXTRMpDPBLD.VRATIO) ... 


» 
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DO 128 1=1,51 
AI«( I-D/50. 

1,23 CAlL DISK VLE , VTE , DPBLD , A I , VRELS ( I ) , VRELP ( I ) , DUDXS < I ) , DUDXP (Dr 
lEXTRM, VRATIO) 

305 DO 306 1=1,51 

PMP2S =PMPIP( I )*( (1 . /TIPDP*( 1 .0-GAmM1*(VR£LS( I>**2) ) > 

1**( l . /GAMM2 ) ) 

PMP2P =PMPI P( I)#((1./TIPDP*(1. 0-GAMMl*( VRELP ( I) ■»*2 ) ) ) 

1 **( 1 ./GAHM2) ) 

IF ( I8LD . EG. 2 ) GO TO 310 
DELP=PMP2P-PMP2S 
GO TO 320 

310 DELP=PMP2S-PNP2P 
I 320 GRAN(I)=HMSTR(I)*RMSTR(I)#DELP 

306 CONTINUE 
GO TO 130 

121 SIGMA=SIG ( 2 ) 

DPBLD “DPBLD+DELBD 
GO TO 130 
125 SIGMA=SIG(2) 

1 WRITE (6, 2129) 

130 CONTINUE 

IF ( DSTRES . NE . 0 . 0 . AND . IBLD . EG . 2 ) CAlL SDISK ( RH3 , RTS , CHORD , RHODI , 

1 RHGBL , RHOAT , SPEED , S I GMA r AREAR , STRESD , STRESB ) 

C ASSUME THE ROTOR SEES THE ROTOR INLET RELATIVE TOTAL TEMPERATURE 
IF(DSTRES.NE.O.O.AND. IBLD- EG, 2) CAlL AL I FE ( STRESB , T2 P PH, SL I FE ) 

* IF( IPRINT.EG. 1 ) WRITE(6,2112> 

IF ( I PRINT. EG . 1 ) WRITE (6,2113) 

IF< IPRINT.EG.O)GO TO 502 
DO 2000 j=l ,51 ,5 

WRITE (6, 21 14) AM( J ) , VRELS ( J) ,VRELP( J) 

2000 CONTINUE' 

> 502 RED IS =1 .-(VRELS ( 1 ) /VRELS (51 ) ) **2 

DSDIS=DPBLD 
DO 131 1=1,51 
131 VDP(I)=VRELP(I) 

CALL SMAlST ( VDP , 51 , MDP ) 

DPDIS=1.-VDP(MDP)/VDP(1 ) 

1 IF ( I PRINT. EG . 1 ) WRITE ( 6 ,2140 > RED I S , DPDIS , DSDI S 

IF ( DSTRES. NE. 0.0. AND. RBEST .EG. 0.0) GO TO 501 
IF ( ITRIG . EG . 0) GO TO 501 

I F ( KZWFL . EG. 1 .AND. IPRINT.EG. 1 ) WRI TE ( 6 , 21 43 ) XWEIFL 
501 TMT2 ( 1 ) = 1 . 0 

VUCRM(1 )=VUCRI 

) WUCRM ( 1 ) = ( VUCRI-RPM ) #SGRT( T I PDP > 

BETAfM 1 ) =BETLE 

EPSM( 1 ) = 1 . O-TLE/ ( SPAC*COS ( BETLE ) ) 

THIK ( 1 ) =RMSTR ( 1 ) * ( 1 . O-EPSM ( 1 ) ) /SIGMA 
ANGCK=TAN(BETTE) 

BNGCK=TAN ( BETLE) 

v WXCR I = VRELS ( 1 ) *COS ( ANGZ 1 ) 

WXCRO= VRELS ( 51 ) #CQS ( ANGZ2 ) 

DO 40 K=2, 51 
KNT =K-1 
AKNT=.(NT 

IF ( KNT. NE. 1 ) GO TO 43 
SUMM=0 . 5* ( GRAN ( 1 ) +GRAN ( 2 ) )*0.02 
1 GO TO 41 

43 SUMM=SIMP2(GRAN, 0.0, 0.02, KNT) 

41 IF ( IBLD .EG .2) GO TO 44 

TMT2(K ) =1 .0 
GO TO 42 

44 TMT2(K)=1 . 0+GAMM2*fiPM*SIGMA#SUMM/ <RHOV( GAMMA , VCRI )#CGS(AIN) ) 

I 42 VUCRM(K)=(VUCRI+(GAMMA+1 . 0) *SIGMA*SUMM/ (2 . *GAMMA*RHOV ( GAMMA , VCRI ) 

1 *COS (AIN) ) ) / ( RMSTR ( K ) *SGRT ( TMT2 ( K ) ) ) 

IF (IBLD. EG. 2) GO TO 45 ... 


TPTPM ( K ) = 1 . 0 

GO TO 46 92 

45 TPTPM ( K ) 3 1 . O-GAmM 1 # < RPM**2 ) /TMT2 ( K ) 4 C RMS7R ( K ) **2 ) * < 2 . *VUCRM { K ) * 
1SGR7 ( TMT2 < K ) ) / ( RPM*RMS7R ( K ) ) -1 . 0 ) 

46 wUCRM <K ) = ( VUCRM (K ) -RPM*RMSTR ( K ) *SGR7 ( i . /TMT2 (K ) ) 1 / SORT < TPTPMiK ) ) 
WXCRM=NXCR 1 + f NXCRO-NXCR I ) *AKNT/50 . 

ARGCK =NUCRM { K ) /WXCRM 

IF( IBLD.EG.2) GO TO 51 

IF < ARGCK .LE.AnGCK) GO 70 52 

BETArt ( K ) a 8E7TE 

WCRM<K >=NUCRM(K> /SIN (BETTE) 

GO TO 53 

51 IF ( NUCRM ( K ) . LT ■ 0 .0 ) GO TG 54 
IF< ARGCK . l £ • SNGCK ) GO TO 52 
BETAri ( K ) =BE7LE 

NCRM ( K ) =NUCRM ( K ) /S IN < BETLE ) 

GO TO 52 

54 IF < A8S ( ARGCK > .LT. ABS< ANGCK > ) GO TO 52 
BE7AM(K>=BE77E 
WORM ( K ) 3 NUCRM ( K ) /S I N < BETTE ) 

GO TO 53 

52 BETAM(K ) =ATAN< ARGCK ) 

NCRM(K) =SGRT<WXCRM**2+NIjC ' )*#2) 

53 EPSM<K>=RHGV(GAMMA,VCRI )*COS(AIN)*PIPDP*SGRT( 1 ./TIPDP>*< 1 ./ 

1PMPIP(K > ) *SGRT (TMT2P ( K ) )/( RhOV( GAMMA , wCRM < K ) ) *HMS7R ( K ) *RMSTR (K ) * 
2C0S ( 8E7AM { K ) ) ) 

THIK (K >=RMS7R(K )*( 1 .0-EPSw(K >> /SIGMA 
40 CONTINUE 

SUM=0. 

CXLEN(1)«1.E-10 
DO 25 1=2,51 

BET AV= f BET AM ( I ) +8ET AM ( I - 1 ) >/2. 

SUM=SUM+ABS (0 . 02/C0S ( BE7AV ) ) 

CXLEN(I)=SUM 
35 CONTINUE 

THIKR( 1 ) =7HIK ( 1 ) /2 . 

THIKL (l)=-THIK(l)/2. 

THETA ( 1 >=0.0 
RTHET ( 1 ) =0 . 0 
R( 1 ) =RIN 

AMAVE=GAMM1*(VELBL**2> 

CALL HM(AMAUErHFAVE) 

ANGMN = ATAN( (TAN<ANGZ1 > +TAN < A«GZ2 > > /2 . ) 

CLSC=2.*<TAN<ANGZl )-TAN(ANGZ2) ) *COS ( ANGMN ) 

SEC0=0.0l38 

YSEC= ( COS ( ANGZ2 ) /COS ( ANGZ 1 ) > * ( 0 . 0075* { CLSC**2 ) * < COS < ANGZ2 ) **2 ) / 

1 < COS < ANGMN ) **3 ) +0 . 035 ) 

PSPTTE= ( 1 . -AMAVE > ** ( 1 . /<3AMM2 ) 

DO 90 1=1,51 

AKS=(l.-GAMMl*< ORELS < I)**2> ) **0.467 
AK P= ( 1 . -GAMM 1 * ( VRELP ( I ) **2 ) ) **0 . 467 
AMS=GAMMl* < ORELS ( I ) **2 > 

AMP=GAMMl * < VRELP ( I ) **2 ) 

CAlL HM(AMS,HMS< I ) ) 
call HM< AMP ,HMP ( I ) ) 

TVIS=TMT2P< I)*TIVIS 
CALL UISCO (TV IS, VIS) 

ASON= ACR I T ( GAMMA , RG AS , TV I S ) 

NS ( I ) = VRELS ( I ) *A30N 
WP ( I ) =URELP < I ) *ASGN 
ZS=NS < I ) 

ZP=NP( I) 

RHQS=RHOV ( GAMMA , ORELS ( I ) ) /MRELS ( I ) *PMP I P < I ) /TMT2P ( I > *RH02 
RHOP=RHOM ( GAMMA , VRELP ( I ) ) /URELP ( I ) *PMPI p ( I ) /TMT2P < I ) *RH02 
ROS ( I ) =RHOS 

ROP(I)=RHOP 


I 

ft 

ft 

ft 

90 

93 

ft 

ft 

ft 

ft 

ft 96 

ft 

» 

OR 

t w 


CXAC7=CHGRD*CXLEN< I ) 

R£0S=RH0S*ZS*CXAC7/VI 5/ 1 2 . 

R£0P*RhGP*ZP*CXAC7/VIS/l2. 93 

EXPOS* 1 . /<2.6*<RE0S**< 1 ./14. ) > ) 

£XPQP°1 . /(2.6*<R£0P**U ./l4. ) ) ) 

IF(R£OS . L7 . 4200. >EXPGS=0.2l2 
I F i REOP . I_7 , 4200 . ) EXP0P=0 .212 
VEX0S=1 ./EXPOS 
VEXOP-1. /EXPOP 

FS(I)«( ( ( ( RHOV ( GAnMA r ORELS ( I ) ) * < ORELS ( I > ** ( I . +HMS (I ) ) ) )**1 .260) 
1*(<VIS/(RH0S*ZS) )**O.26S)*AKS)/(l0.**(0.678*(2.*EXP0S+l . ) ) )) 

2/COS <B£7AM< I ) ) 

FP (I ) o t ( ( < RHOV ( GAMMA , VRELP < I > )■*< VRELP ( I > **< 1 . +HMP ( I >)>)*# 1 . 268 ) 

X*< (VIS/ ( RHOPfl'ZP) >**0.263)*A«P>/< 10 . ** < 0 . 678* < 2 . *£XP0P+1 . ) ) ) ) 

2/C0S (BE7AM( I > ) 

FSL <!)=<( VRELS ( I > /VELBL ) **5 ) /CDS ( BETAn (I) ) 

FPL ( I ) = ( ( VRELP ( I ) /VELBL ) **5 ) /COS ( BE7AM ( I ) ) 

BNUS ( I ) *VIS/RHQS 
BNUP( I ) »VIS/RH0P 

DUOXS < I ) = l 2 . *DUDXS ( I ) *AS0 N*CGS ( BETAm ( I > ) 

DUDXP ( I ) = 12 . *DUDXP < I ) *ASON*COS ( BETAM ( I ) ) 

C0N7INUE 

92 CHRD ( ICHORD ) =CHGRD 

REYN*RHORE*VELRE*SUM*CHORD/ ( 1 2 . *V I SRE ) 

RNGSE*7LE/24. 

REZNO*NS ( 1 ) *RNQSE/BNUS ( 1 ) 

IF < RNOSE . NE . 0 . 0 ) 7H£S0=0 . 21 *RN0SE*1 2 . /SGR7 ( REZNO ) 

I F ( RNOSE . EG . 0 . 0 ) 7H£S0=0.0 
7HEP0 3 7HESO 
CAML7=CH0RD#SIJM 
XL£N*0 . 02 

CFLAM= 1 . G28/SGR7 ( REYN ) 

INGS = 1 
INGP= 1 
I FLAGS =0 
I FLAGP=0 
DO 95 1-2/51 
J = I -1 

IF( 1 FLAGS . EG. 1 ) GO 70 96 
INGS=ING5+1 

7HES J = < ( VELSL/VRELS ( I ) > **2 > *0 . 5*CFLAM*SGRT< SIMPZ (FSL ,0 . 0 .XLEN / J ) > * 
1CAML7+7HES0 

ZTRUKS* ( ( 7HESJ/12 . ) **2 ) /BNUS ( I ) 

F7RUKS=Z7RUKS*DUDXS ( I ) /CHORD 
FLAMS S FI7 (F7RUKS) 

HLAnS" (0.3-1 ./l20.*FLAMS)/(37./2l5.-l . /945 . *FLAriS- 1 . /9072 . * 
l(FLArtS**2> ) 

DELSJ=HLArlS*THESJ 

BLREYS=wS ( I ) *DELS J / 1 2 . /BNUS ( I ) 

RECR7S=RECRIT< FLAMS ) 

IF ( BLREYS . GE. RECR7S > I FLAGS* 1 
IF( IFLAGP.EG. 1 ) GO TO 95 
INGP=INGP+1 

THEP J => ( ( VELBL / VRELP < I ) >**3)*0. 5*CFL Ar1*SGR7 < S I MP2 ( FPL , 0 . 0 , XLEN , J > ) * 
1CAML7+7HEP0 

ZTRUKP=( (7HEPJ/12. > **2 > /BNUP ( I > 

F7RUKP*Z7RUKP*DUDXP< I ) /CHORD 
FLAMP=FI7 (F7RUKP) 

HLArtP" ( 0 . 3-1 . / 1 20 . *FLAnP > / ( 27 . /31 5 . -1 . /945 . *FLAMP-1 . /9072 * 

1 <FLAMP**2> > 

DELP J «HLAftP*7HEP J 

BLR£YP=WP< I )*DELPJ/12./BNUP< I ) 

RECRTP=RECR I T ( FLAMP ) 

I F < BLREYP . GE . RECR7P ) I FL AGP* 1 
CGN7INUE 

YL£N=CH0RD/600. 
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SPACE=CHQRD/SIGMA.*CGS(B£7Ar!(5l > ) 

NS=52-INGS 
NP=S2-INGP 
IS=NS-1 
IP=NP-1 
DO 97 1=1 ,NS 
J=INGS*I-1 
GS( I ) =FS < J ) 

97 CONTINUE 
DO 96 1=1 ,NP 

j J= INGP+I-1 

GP( 1 ) =FP( J ) 

98 CONTINUE 

TH£S2=< (0. l56/( (RH00(GArtMA,VR£LS(51 ) )*(URELS(51 >**< 1 .+HMS(5l ) ) ) ) 

J 1**1 .268) *SIMP2(GS, 0.0, YLEN, IS ) + ( ROS ( INGS ) * (NS ( INGS ) ** ( 2 . +HMS ( I NGS ) 

{ 2) )/(R0S(51 >*(NS(51 )**(2.+HMS<51 ) ) ) )*(TH£SJ/l2. ) )**1 . 266) **0. 7386) 

2 12. 

THEP3=( <0. l56/( ( RHOV( GAMMA, VR£LP( 51 ) ) * ( ORELP ( 5l >**( 1 .+HMP(5l ) > ) ) 
1**1 .268>*SIMP2(GP,0.0, YLEN, IP)+<ROP( INGP)*(HP( INGP > ** ( 2 . +HMP ( INGP) 
;i 2) )/<R0P(51 )*<HP(51 >**(2.+HMP(5l ) ) ) )*<THEPJ/1Z. ) )** i . 268 ) **0 . 7866 ) 

3*12. 

TSTAR= ( THES3+THEP3 ) /SPACE 
DELS3=HFAVE*THES3 

;i IF ( I FLAGS . EG . 0 ) DELS3=DELS J 

D£LP3=HFAME*THEP3 

■ j IF< IFLAGP . EG . 0) DELP3=D£LP J 

) DSTAR= ( DELS3+DELP3 ) /SPACE 

DTE{ 1 ) =TETS/SPAC£ 

DTE (2) =TETR/SPACE 
AFSl°GA(1Ml*(NCRM(5l )**2) 

BLOK 1 = 1. -DST AR-DTE ( I BLD ) -TST AR 
BLQK2=1 .-DSTAR-DTE(IBLD) 

CSTN= ( ( 1 . -AFSl )*( (GAMMA+1 . ) / < 2 . *GAMMA ) )+8L0Ki*( (NORM (Si )* 
lC0S(B£TAM(51 ) ) )**2) ) / (BL0K2*fcCRM ( 51 ) *CGS (BETAM < 51 ) ) ) 

DSTW=BLQK 1*HCRM( 51 ) *SIN ( BET Ad (51 ) )/BLGK2 

WXM IX=GAMMA*CSTW/ (GAmMA+ 1 . )-SGRT( < ( GAMMA*CSTW/ ( GAMMA* 1 . ) )**2)-l .0+ 
1 GAMM 1 * ( DSTW**2 ) ) 

RHMIX=( 1 .-GAMM1*(HXMIX**2+DSTW**2) )**GAMM3 

■ PLOSS ( IBLD) =BL0K2*RHGV (GAMMA, wCRM( 51 ) ) *COS ( BETAM ( 51 ) >/(RHMIX* 
1WXMIX) 

SPLOSS=PLOSS( IBLD) 

YPRO ( IBLD ) = ( 1 . -PLOSS ( IBLD ) ) / ( PLOSS < I BLD ) * < 1 . -PS PTTE ) ) 

YSF ( IBLD) =YSEC*CH0RD/HEIGH7 
YTOT ( IBLD ) = YPRO ( IBLD ) +YSF ( I BLD ) 

PLOSS ( IBLD ) = 1 . / ( 1 . +YTOT ( IBLD ) * ( 1 . -PSPTTE) ) 

IF< ITRIG.EQ. 1 )GQ TO 94 
YTOTL( I CHORD ) =YTDT ( IBLD ) 

I F ( I CHORD . LE . 2 ) GO TO 400 

I F ( YTOTL ( I CHORD ) . GT . YTGTL ( I CHORD- 1 ) . AND . YTOTL ( I CHORD- 1 ) .07. 

1 YTOTL ( ICHQRD-2) ) GO TO 91 
400 I CHORD® I CHORD* 1 

CHGRD=CH0RD+0 . 02 
IF( CHORD. GT.CHMAX) GO TO 91 
IF(TTE.tG.O.O) GO TO 92 
GO TO 500 

91 MI NMUM= I CHORD-2 

CHORD =CHRD (MI NMUM ) 

’ ITRIG= 1 

KNZN= 1 

IF(TTE.EG.O.O) GO TO 92 
GO TO 500 
94 CONTINUE 

TRANS =CXL£N ( I NGS ) *CHORD/CAMLT* 1 00 . 

: ) TRANP=CXLEN ( I NGP ) *CHGRD/CAMLT * 1 00 . 

ALPHA=RIN*(R3STR-1 .) /CHORD 

! . _D.O .70 N=2.,S1 : 
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NN=N-1 

T0= THETA (NN> 

R<N)=RIN*RMSTR<N> 

BE=BETAM(NN> 

IF<BE.EG.0.0> BE=0. OOOOOl 

76 RHS 1 =R < NN ) *THETA < NN ) + < TAlM ( BE ) +TO*Al PHA > *0 . 02*CH0RD 

75 700=RHSl/R<N> 

RHS2 = R < NN > *THE7A ( NN > + < TAlM ( BE ) +T00 *AlPHA ) *0 . 02*CH0RD 
IF ( ABS ( (RHS2-RHS1 >/RHSl >-.00002) 71 ,71,72 

72 KN=KN+1 
IF (KN-100) 73,73,74 

73 RHSl =RHS2 
GO TO 75 

74 IF< IPRINT . EG. 1 > WRITE(6,2l IS) 

71 THETA (N)=RHS2/R(N> 

RTHET (N) =R (N ) *THETA < N ) /CHORD 
THIKR(N) » RTHET <N)+THIK <N) /2. 

THIKL(N)=RTH£T(N)-THIK(N>/2. 

70 CONTINUE 

PITCH=1 . /SIGMA 
D I STAN 01 ABS ( RTHET <51 > > 

STAGR=ATAn(DISTAN> 

STAGR=STAGR* 1 80 . /3 . 1 4 1 5? 

CAMLTH=CHGRD»SUM 
ZNUM=2.»3.l4lS9*RIN*SIGMA/CH0RD 
NUMBER=ZNUM 

IF ( I PRINT . EG . 1 ) WRITE (6, 2115) 

IF < I PR I NT. EG. 1 ) WRITE <6,2116) 

TCAL=THIK <51 >*COS (BETAM'51 ) )*CHORD 
ANACT=BETTE*l80./3. 14159 

101 T3IT=TMT2(51 )*T1P 
ANIT=BETAM<51 >#lS0./3. I4i59 
IF( IBLD.EG.2) GO TO 102 
IF ( I PRINT. EG. 1) WRITE ( 6 » 21 17 ) T2IT ,T1P, UUCRM ( 5 1 ) , VUCRO , AnIT, ANACT, 

1 SIGMA, SOLID, TCAL,TETS 
GO TO 103 

102 IF < IPRINT. EG. 1 > WRITE (6 , 21 17 ) T3IT , T3P , VUCRM<51 ) ,UUCRG,AnIT, AnACT , 
1SIGMA, SOLID, TCAL,TETR 

103 CONTINUE 

IF ( IPRINT .EG . 1 ) WRITE ( 6 ,2100 ) 

IF< IBLD . EG . 2 ) GO TO SO 

IF ( I PRINT . EG. 1 ) WRITE ( 6 , 21 1 9 ) 

GO TO SI 

80 IF( I PRINT . EG . 1 ) WRITE ( 6 , 21 20 ) 

81 IF< IPRINT .EG. 1 ) WRITE(6,2121 ) 

DO 82 1=1,51,5 

IFC IPRINT. EG. 1 ) WRITE (6, 2 122) AM ( I ) ,R(Ii ,RTH£T(I) ,THIKR(I) , TH I K L < I ) 

82 CONTINUE 

IF ( IPRINT .EG . 1 ) WRITE (6 ,2123 ) 

I F ( I PR INT . EG . 1 > f ,R I TE < 6 , 2 1 24 ) CHORD , S I GMA , CAMLTH , STAGR , NUMBER 

IF( IBLD .EG .2 ) GO TO 110 

IF< I PR INT. EG. 1 ) WRITE(6,2125) 

GO TO 111 

110 IF( IPRINT. EG. 1) WRIT-' (6, 2 126) 

111 I F < I PR INT . EG . 1 ) WRITE(6,2127) 

UMIX°SGRT ( WXMIX##2+DSTW*#2 ) 

DSTW=A8S( DSTW) 

IF( IPRINT. EG. 1) WRITE ( 6 , 212S ) TSTAR ,DSTAR ,OMIX , SPLOSS 

IF( IPRINT. EG. 1) WRITE (6, 2144 )YPRO( IBLD) ,YSF< IBLD) ,YTOT< IBLD) , REYN 

IF( IPRINT. EG. 1 ) WRITE(6 ,2221 ) TRANS, TRAnP 

IF (IPRINT. EG. 1 ) WRITE (6, 2 100) 

IBLD=IBLD+1 

IF( IBLD. LT . 3 ) GO TO 20 

ETATT = 1 . / ( 1 . 0+ ( T3P/ ( T1 P-T3P > * < 1 . 0- < < PLOSS ( 1 ) *PLOSS < 2 ) ) **GAMM2 ) ) ) > 
IF( IPRINT. EG. 1) WRI7£(6 , 2099) .... _ 



IF( I PR I NT . EG . 1 ) WRITE ( 6 f 2122 ) 

IF< IPRINT. EG . 1 ) WR±T5(6f2133) ETATT 

VX2C=VX3*ACR3 96 

VX2C=V2*ACR1 *CGS ( AlP2 ) 

UUrt2C=VWRK3 

UM3C=U2 

VUM2C=UWRK2 

UM2C°U2 

R3C = i .+ h3/(2.*RM3> 

R2C a 1 . +H2/ < 2 . *RM2 > 

R3CI =1 . /R3C 
R2CI=1 ./R2C 

REACTC = ( UX3C**2+ ( R3CI *UUM3C-R2C*UM3C > **2-UX2C**2- < R2CI *VUN2C-R2C 
1 *UM2C ) **2 ) / < VX3C**2+ ( RSC I *UUH3C-R2C*UN3C ) *#2+2 . *Un2C*VUM2C- ( R2C 
2*UM2C>**2) 

HROTOR® ( H2+H3 ) *6 . 

ETA7TM= ETATT* ( 1 . -CLEAR /HROTOR* ( 2 . 755*( R£ACTC**2 5 +0 . iOS*REACTC+ 
11.72)) 

DELETA»ETATTh-ETATT 

IF< IPRIMT.EG.l) WRITE ( 6 f 2l45) DELETE . tTATTM 
PTTEX«P1P*( ( 1 .-( 1.-T3P/T1P) /ETATT) **( 1 . /GANM2) ) 

PSTEX=PTTEX* ( ( 1 . -GAUM l * ( 03**2 ) ) ** ( 1 . /GAnM2 ) ) 

PXTEX = PSTEX/ ( ( 1 . -GAllft 1 * ( UX3**2 ) ) ** < 1 . /GAilM2 ) ) 

ETATX= < ETA* ( 1 . - ( PTTEX/Pl P ) *#GAf1fl2 ) / ( 1 . - ( PXTEX/Pl P ) **GAMN2 ) ) +DELETA 
ETATST=(ETA*< 1 .-( PTTEX/Pl P ) **GAflN2 ) / ( 1 . - < PSTEX/Pl P > **GAflM2 ) ) + 
1DELE7A 

IF( IPRINT .EG. 1 ) WRITE ( 6 f 2220 ) ETA VET ,ETATX 
IF ( DSTRES. NE. 0.0) GO TO 600 
HUBRa ( RM2-H2/2 • )*12. 

TIPR* (RM2+H3/2 . )*12. 

REARA a 3 .141 59* ( T I PR**2-HUBR**2 ) 

CALL SDISK (HUBR fTIPR fCHORD fRHODI fRHOBL fRHQAT f SPEED fS IGMA f REAR A f 
ISTRESDfSTRESB) 

IF(I PRINT .EG . 1 ) WRITE ( 6 f 224 1 ) ETREED , STRESB 
600 SEFFY ( NPASS ) a E7ATST 

PCHEK =A8S ( ETATT-ETA ) 

I F ( PCHEK . LE . 0 . 002 . AND . DSTRES . EG . 0 . 0 ) I PR I NT= 1 
IF ( PCHEK . J_E . 0 . 002 . AND . RBEST . GT . 0 . 0 ) I PR I NT = 1 
IF (NPASS. EG. 1 ) GO TO 150 
DETA=ASS ( SEFFY (NPASS ) -SEFFY ( NPASS-1 ) ) 

IFtDETA.LE. 0.002) GO TO 142 
IF ( NPASS . GT . 6 ) GOTO 141 
150 NPASSaNPASS+1 
E7A=ETATT 
P2P=PLQSS(1 )*P1P 
PLOSl = PLOSS ( 1 ) 

PL0S2 a PL0SS(2)*(T22PP**(l ./GANM2) ) 

P32PP=PL0S2 

IF ( DSTRES. NE .0.0) GO TO 180 
GO TO 10 

141 WRITE(6f2124) 

142 CONTINUE 

IF] ( RBEST . GT . 0 . 0 ) WR I TE < 6 f 2234 ) DSTRES f STRESD f STRESB 

IF (IPRINT. EG. 1)WRITE(6, 2142) NPASS 

IF (DSTRES .EQ.O .0) GO TO 4000 

I F ( RBEST . GT . 0 . 0 ) GO TO 4000 

KPASS=KPASS+1 

EADIAB(KPASS) =1 . -E7ATTM 

DPSI (KPASS) =STRESD 

UHUBC a UT3*RH2/RT3 

I F ( K PASS . EG . 1 ) WR I TE ( 6 f 2 1 00 ) 

IF(KPASS.£G. 1 ) WR ITE ( 6 f 2230 ) 

IF (KPASS . EG . 1 ) WRITE i 6 f 2244 ) 

IF(KPASS.EG. 1 ) WRITE (6f 2231 ) 

WR I TE ( 6 , 2232 ) K PASS f E7A7TM f ETATST f UHUBC , UTS , RH3 , RT3 f STRESB f STRESD , 
1BLIFE 


IF(V2.G£. 1 .0>WRITE(6*2237> KPASS 
IF(WPP3.GE.l .0>WRITE(6*222S) KPASS 

IF (KPASS . L£. 20) GO TO 1?0 97 

K =0 

DO 450 1=1*21 

IFfDPSI < I > .GT.DSTRES) GO TO 450 
K = K + 1 

FADIA8(K > =EADIAS( I ) 

RADEXX < K > = R ADEX { I ) 

CONTINUE 

IF(K.EG.O) GO TO 460 
WRITE (6, 2236) K 
I F ( K . EG . 1 ) nlNEf A=K 

IF(K.GT. l )CALL SMAlST (FADIAB.K *MlN£TA) 

RB£ST=RAD£XX ( MINETA ) 

KPASS=0 
GO TO 190 
WRITE (6 * 2235 ) 

GO TO 4000 
FORMAT ( 7F 10.3) 
t FORMAT ( 6F10 . 3 * 1 1 ) 

2 FORMAT (3F1 0.3) 

7 FORMAT < /// ) 

FORMAT < ///56X * 1 1HFL0W ANGLES/56X , 1 1H i 

FORMAT <////////) 

FORMAT ( 1H1 ) 

FORMAT ( SOX * Z4HPR0GRAM INPUT PARAMETERS/50X , 24H 

1 ///> 

FORMAT ( SX* 5HGAMMA* 6X * 4HRGAS * 6X * 11HRGT0R SPEED * 6X * 1 4HMASS FLOW RATE 
1 *6X* l3HTUReiN£ INLET *6X* 13HTUR8INE INLET * 6X * 1 3HW0RK RE5UIRED/32X * 
25H(RPM) * 12X,SH (LB-SEC) ,9X* 14HPRESSURE (PSI ) *5X* 1 4HTEMPERATURE * R, 
38X * SH ( STU/LS ) / ' 

FORMAT £ 7X * F6 . 3 * 4X * F7 • 3 * 7X * F7 . 1 , 1 lX * F7‘. 3 * 13X * F7 . 3 , l2X ,F7 . l , 12X*F7.3 
1 ) 

FORMAT (///////2X* 12HSTAT0R INLET ,2X * 1 1HR0T0R INLET *2X * 1 OHROTOR EXI 
lT,2X* 12HSTAT0R INLET * 2X * 1 1HR0T0R INLET , 2X , 10HR0TGR £XIT*2X* 1 1HSTAT 
20R T , £ . * 2X* 1 OHROTOR T . E . * 2X * 1 5HR0TGR CLEARAnCE/ZX * 1 1 kMEAN RADIUS* 
33X * 1 1 HME AW RAD I US * 2X * 1 1 HflEAN RAD I US * 5X , 4H AREA , 1 OX * 4H AREA , SX * 4HAREA 
4 * 6X * 9HTH I CK NESS , 3X * 9HTH I CK NESS * 5X * SH ( I NCHE5 ) /4X * SH ( I NCHES ) , 6X * SH ( I 
5NCHES) *4X,SH< INCHES) *3X, llh<SG INCHES) *3X, 11H(SG INCHES ) , 2X , 11H ( SG 
6 INCHES) *3X,SH( INCHES) , 4X, SH ( INCHES )/ ) 

FORMAT < 5X * F6 . 3 , SX * F6 . 3 * 6X * F6 . 3 * 7X * F7 . 3 * 6X * F7 . 3 * 6X * F7 . 3 * 6X * F6 . 3 * 6X * 
:F6.3,9X*F6.3) 

FORMAT < //// 40X * 46HC0MPUTED FLOW PARAMETERS AND VELOCITY DIAGRAMS/ 

140X * 46H ///) 

FORMAT • 22X * 6HSTATGR * 60X * 5HR0T0R/22X , 6H , 60X , 5H / / /2X * 

1S0HSTAT0R INLET CRITICAL VELOCITY RATIO* V/VCR ,F7.3*5X* 

ZSOHRQTOR EXIT AXIAl CRITICAL VELOCITY RATIO* VX/VCR . ,F 7.3//2X, 

350HSTAT0R INLET SWIRL VELOCITY RATIO* VU/VCR . ,*F7.3*5X* 

450HR0T0R EXIT SWIRL VELOCITY RATIO* VU/VCR ,F7.3//2X* 

550HSTAT0R EXIT CRITICAL VELOCITY RATIO* V/VCR ....... *F7. 3, 5X* 

650HR0T0R EXIT CRITICAL VELOCITY RATIO, V/VCR ,F7.3//2X, 

750HSTAT0R EXIT AXIAL VELOCITY RATIO* VX/VCR ,F7.3*5X* 

350HR0T0R INLET RELATIVE VELOCITY RATIO, w/WCR *F7.3/) 

FORMAT (2X*50HSTAT0R EXIT SWIRL VELOCITY RATIO* VU/VCR * 

1 F7 . 3 * 5X * 50HRQT0R EXIT RELATIVE VELOCITY RATIO, w/WCR *F7.3 

2//2X *49HSTAT0R EXIT ABSOLUTE TOTAL PRESSURE *F8.3* 


F7 . 3//2X , 

...... *i-7. 3* 5X * 

...... *1-7. £/ / 2X * 

......*i s 7.w* 5 X * 

. . . . . . *F7.3/) 

VU/VCR * 

w/WCR * F7 . 3 

»FS .3* 


35X , 50HRQTQR INLET BLADE SPEcD RATIO* U/VCR .*F7.3//2X* 

450HSTAT0R ABSOLUTE TOTAl PRESSURE LOSS RATIO* ,F7.3*5X* 

550HR0T0R EXIT BLADE SPEED RATIO* U/VCR *F7.3//64X* 

650HR0T0R RELATIVE TOTAl PRESSURE LOSS RATIO *F7.3//64X* 

749HR0T0R EXIT ABSOLUTE TOTAL PRESSURE *F8.3//64X* 

349HR0T0R EXIT ABSOLUTE TOTAL TEMPERATURE ,FS.1//64X* 

949HRGT0R TIP SPEED (FT PER SEC) *F8.1> 

FORMAT (//59X*5HSTAGE/59X*5H ///37X*45HR0T0R EXIT HUB RADIUS (I 

INCHES) *F7.3//37X*45HR0T0R EXIT TIP RADIUS (INCHES) 



2 fF7.2//27Xf45HR0TGR INLET HUB RADIUS (INCHES) ...... 

2 7 . 3//37X » 45HRGT0R INLET TIP RADIUS (INCHES) 

4 » F7. 3//37X » 45HSTA7GR INLET HUB RADIUS (INCHES) ,F7.3// 98 

537X, 45HSTAT0R INLET TIP RADIUS (INCHES) . ,F 7.2//37X, 

645HR0T0R EXIT ANNULUS AREA (SG INCHES) , F7. 3//37X . 4SHR0T0 

7R INLET ANNULUS AREA (Sl3 INCHES) F7.3//37X.45HS7AGE REACT 

SIGN F 7.3/) 

2110 FORMAT (////36X fSZHCAlCUL ATIONS FOR GENERATION GF STATOR BLADE GEGn 

1 ETRY/36X . 52H > 

2111 F0RMAT(////36X,51HCAlCULATI0NS FOR GENERATION OF ROTOR SLADE GEGME 

l TRY/36X . SI H > 

2112 FORMAT (//44X,36H**» COMPUTED AERODYNAMIC LOADING »#»/> 

2113 FORMAT ( 30X » ' PERCENT SUCTION SURFACE 

1 PRESSURE SURFACE '/29X # ' MERIDIONAL RELATIVE CRITICAL 

2 RELATIVE CRITICAL '/30X, ' DISTANCE VELOCIT 

3Y RATIO VELOCITY RATIO'/ 5 

2114 F0RMAT(3iX,FS.3f21XfF5.3,23XfF5.3) 

2115 FORMAT < ////3SX,55H*** (ITERATED) STAnITZ METHOD SLADE EXIT PARAMET 
1ERS **#/) 

2116 FORMAT (45X , 1 4HITERATED VALUE .SX f IZHACTUAl VAlUE/) 

2117 FORMAT ( 1 1 X r 27HEX I T (ABSOLUTE) TEMPERATURE , 9X, F9. 4 r l2X,F? . 4//3X, 
135HEXIT (ABSOLUTE) TANGENTIAL VELOCITY f 1 IX ,F7 .4 f 14X fF7 . 4//22X , 

Z16HEXIT BLADE ANGLE f 1 OX f FS . 4 f 1 3X f FB . 4//30X , gHSOLIDITY f 1 1 X r 

3F7 .4 f 14X f F7 . 4// 15X. 23HTRAILING EDGE THICKNESS , 12X ,F6 .3 , 1SX rfi . 3) 

2118 FORMAT ( /7lH **** RflTHETA ITERATION LOOP NOT SATISFIED CQMPUTA7I0 

IN CONTINUES ) 

2119 FORMAT ( / ///22X , 4H*#-» f 12Xf 26HFINAL STATOR BLALt PROFILE, 13X,3H*»*/ 
133X,S8H*** COORDINATES NORMALIZED niTH RESPECT TO BLADE CHORD ***/ 

2) 

2120 F 0RMAT(////33X,4H*** f 13Xf25HFINAL ROTOR BLADE PR0FILEf12X,3H***/ 
133Xf58H**« COORDINATES NORMALIZED WITH RESPECT TO SlADE CHORD ***/ 

2) 

2121 FORMAT ( 3QX f 7HPERCENT , lOX f 9HP ITCHLINE f 10X f 4HMEAn f 10X f 5H8LADE f IOXf 

1 5HBLADE/29X f 1 OHMERIDIONAl f 9X f 6HRADIUS , 9X , 1 OHCAMBERL INE f 6X f 7HSURFAC 
2E f 8X , 7HSURFACE/30X f SHD I STANCE f 25X , 1 OHCOORDINATE f 5X f 1 OhCOORD I NATE f 
35X , 1 OHCOORD I NATE/ ) 

2122 F0RMAT(31XfF5.3f1ZXfF6.3f10XfF7.4f8X,F7.4f3XfF7.4) 

2123 FORMAT ( //21Xf llhBLADE CHORD ,6X, 14HBLADE SOL I D I TYf 6X f 1 6hBLAD£ CAmBE 
1 RL INE f 6X f 1 3HBLADE STAGGER f 6X , 1 2HBLADE NUMBER/22X f SH ( I NCHES ) f 2SX f 1 5 
^HLENGTH ( INCHES) » 

21 IXfSH ANGLE/) 

2124 F0RMAT(23XfF7.4,UX»F8.4,i3XfF8.4,i3XfF6.2f14X, 13) 

2125 FORMAT(/// 23Xf' CALCULATIONS FOR GENERATION OF STATOR SLADE RON 

1 AERODYNAMIC PERFORMANCE LOSSES ' /23X f ' 

2 ' j 

2126 FORMAT ( /// 23Xf ' CALCULATIONS FOR GENERATION OF RGTOR BLADE RON A 

lERQDYNAMlC PERFORMANCE LOSSES '/22X, ' 

2 F } 

2127 FORMAT ( //43X,38H*#* STEWART MIXING LOSS PARAMETERS *»*/ ) 

2128 FORMAT (32Xf40HT0TAl MOMENTUM THICKNESS ( DIMENSIONLESS ) f 1 4X f F 6 . 4/ 
l32XF44hT0TAL DISPLACEMENT THICKNESS (DIMENSIONLESS ) f 10X,F*. 4/ 

232X f 4 1 HCASCADE EXIT (MIXED) CRITICAL MACH NUMBER, 13X,F6, 4/ 

332X , 2SHPR0F I LE (FRICTION) TOTAl PRESSURE LOSS , 16X ,F6 .4 ) 

2129 FORMAT ( ///75H*#-* eLADE LOADING ITERATION LOOP NOT SATISFIED COM 
1PUTATI0N CONTINUES ***///) 

2122 FORMAT ( //// 44X,39HFINAL CALCULATIONS FOR STAGE EFFICIENCY/44X , 

139H > 

2133 FORMAT ' / / /35X , 43HSTAGE ADIABATIC EFFICIENCY CALCULATED FROM BLADE/ 

1 35X , 52HB0UNDARY LAYER AND SECONDARY FLON LOSSES ,F6.4) 

2134 F0RMAT(///3ZX,62H*** NARNING *** STAGE EFFICIENCY ITERATION FAILS 
ID TO CONVERGE) 

2140 FORMAT (//38X,39H0VERAlL BLADE REACTION (R=1-VIN/V0UT> =,5X,F9.3/ 
138X,44HPRESSURE SURFACE DIFFUSION (DP=1-VMIN/VIN) =fF9.3/TSXf44HSU 
2CTI0N SURFACE DIFFUSION (DS«1-V0UT/VMAX> «*»F9.2/> 

2141 FORMAT (//44X ,26HSTAT0R INLET ANGLE ,F9 .3//44X ,26HSTATQR EX 

ITT AWGLF ,FO.3//44X..?AH.RnT0R .Thil.FT. ANG.I.F. . .F9 ., 3 /,/ 


244X,26HR0T0R EXIT ANGl.£ , F 9 . 2//44X , 26HR0T0R EXIT SWIRL AN 

| 2GLE . . . , F9 . 3 > 

2142 FORMAT <////////41X,27H£FFICIENCY ITERATION CONVERGED AFTER ,12, 99 

17H PASSES) 

2143 FORMAT ( /27X , 75HNGTE ... SASIC LOADING DISTRIBUTION COULD NOT SATI 
1SFY SOLIDITY REGUI REMENT/27X » S6hF0R OPTIMUM ZWEIFEL COEFFICIENT .. 

2.. ZWEIFEL CHANGED T0,F5.Z> 

| 2144 FORMAT </46X,33H*** CASCADE LOSS COEFFICIENTS »**//40X , 24HPR0FILE L 

10SS COEFFICIENT, lSX,F6. 4/40X,3lHS£C0NDARY FLOW LOSS COEFFICIENT , EX 
! 2 , F6.4/40X ,3OHT0TAi_ CASCADE LOSS COEFFICIENT ,9X ,F6 .4///29X ,S3HN0TE 

3 CASCADE GEOMETRY OPTIMIZED AT REYNOLDS NUMBER =,£9.2) 

! 2145 FORMAT ( / / 35X , S2HST AGE EFFICIENCY DECREMENT FOR ROTOR CLEARANCE ... 

1 , Fo. 3///35X, 52HFINAC STAGE TOTAL -TOTAL EFFICIENCY 

I 2. ,F6.4) 

j 2146 FORMAT ( 37X , 45HSTAGE WORK COEFFICIENT BASED ON MEAN RADIUS . ,F 7.3// 
127X,45HSTAGE WORK COEFFICIENT BASED ON HUB RADIUS . . ,F7.2//27X‘ 

| 24SHTURBINE TOTAl-TQTAl PRESSURE RATIO ,F 7.2/) 

i 2200 FORMAT <////! OX, 37H** WARNING »# INPUT STATOR EXIT ANGLE IS TOO L 
> lARGE TO SATISFY CONTINUITY REGUIREMENTS//26X, 17HANGLE CHANGED TO , 

\ 2F6 . 2 , 1 X , 7HDEGREES ) 

“ 2201 FORMAT (/////// ?Xrl2HSTAT0R INLET , 6X , 1 1 HSTATOR EXIT, 6X, 14HSTAT0R Z 

j 1WEIFEL , 6X , 13HRGT0R ZWEIFEL, 6X, 12HSTAT0R CHORD , 6X , 1 1HR0T0R CHORD/lO 

I *X, 

5 210HFLDW ANGLE,7X,10HFL0W ANGLE ,3X , 1 1HC0EFFICIENT , 9X , 1 1HC0EFFIC IENT 

• 3, 9X,8H( INCHES) ,9X , 3H ( I NChES ) / ) 

i 2202 FORMAT ( 12X , F6. 2, 12X , F5. 2 , 1 3X, F5. 3 , 15X, F5 . 3 , 14X , F5. 2 , 1 IX, F5 . 3// ) 

1 2203 FORMAT { / / /43X , 22HPR0GRAM DEFAULT VAlUES//21X,54HSTAT0R EXIT FLOW A 

INGLE CALCULATED FROM W0RK/65X , 1 9HAND FLOW CONDITIONS/ 

23lX , 27HSTAT0R ZWEIFEL COEFFICIENT O.S //3l X , 37HR0T0R ZWEIFE 

2L COEFFICIENT 0.8//3lX,S4HSTAT0R CHORD 

4. OPTIMIZED ON MIN IMUM/65X , Z4HCASCADE LOSS C0£FFICIENT//2lX,54HR0T 

50R CHORD OPTIMIZED ON MINIMUM/65X,24HCASCAD£ 

» 6 LOSS COEFFICIENT) 

2204 FORMAT (/////10X,1 1 HSTATOR L.E. , 3X, 10HR0T0R L.E. ,8X,9HAlL0WABLE,3X, 

• 1 1 OHROTOR EX I T , GX , 1 7HBLADE PLOT OPTION , SX, 1 OhROTOR EXIT/1 IX , 9HTHICK 
2NESS , 9X , 9HTHICKNESS , 3X ,11 HDISK STRESS , 7X , 1 OHAXIAl mACH, 12X, 1 1H0 = 

3N0 PLOT, 10X, 11HSWIRL AnGl£/1 2X, 8H < INCHES ), 10X , 2H ( INCHES ), 1 1 X , 5H< PS 
41 ) , 1 2X , 6HNUMBER , 14X , SHI = PLOT// 13X , F5 .3 , 1 2X , F5 .3 , 1 2X , F7 .0 , 1 IX , 

| SF5.3,l9X,Il,l3X,F6.2/> 

j 2210 FORMAT < / /23X , 65HN0TE .... MAXIMUM ALLOWABLE STATOR ANGLE FOR INPUT 
: 1 CONDITIONS IS ,F5.2,3H DEGREES) 

2220 FORMAT ( //3SX , 52HFINAL STAGE TOTAL-STATIC EFFICIENCY 

1 . . ,F6.4///35X,52HFINAl STAGE RATING EFFICIENCY 

Z. . . ,r6 . 4 > 

| 2221 FORMAT < //3X , ' NOTE ..... PREDICTED SUCTION SURFACE BOUNDARY LAYER 

1 TRANSITION AT ' , F5. 1 , * PERCENT OF CAMBER LENGTH '/19X, ' PREDICTED 
2PRESSURE SURFACE BOUNDARY LAYER TRANSITION AT ' , F4. 1 , ' PERCENT OF 
3CAMBER LENGTH') 

2230 FORMAT < // 1 5X , 85HTHE FOLLOWING IS A LIST OF POSSIBLE CONFIGURATIONS 

1 FDR THE INPUT CONDITIONS SPECIFIED/lSX,35H- 

2 //// 

• *) 

2221 FORMAT ( 2X , 4HPASS , 3X , 7HETA T-T,3X,7H£TA T-S,3X,5HU HUB, 3X, SHU TIP, 

12X,1 OHROTOR EXIT , 3X , 1 OhROTGR EXIT, 2X, 12HBLADE STRESS , 3X , 1 1 HDISK ST 

• 2RESS ,3X , lOHROTOR LIFE/29X , 5H ( FPS ) , 2X , 5H < FPS ) ,2X , 1 OhHUB RADIUS, 2X , 1 

*0HTIP RADIUS, 6X, 

■ 35H ( PS I ) ,10X,5H(PSI ) , 7X , 7H< HOURS ) /46X , 3H ( INCHES ) , SX, SH ( INCHES )// ) 

f 2222 FORMA'" X , 12 , 4X ,F6 . 4 , 4X , F6. 4, 4X , FS. 0 , 3X ,F5. 0 , SX, F6. 2, 7X,F6. 2 , SX, 

1F7 . 0 , 7X , F7 . 0 , SX , £ 10 . 4 ) 

2233 F0RMAT(////SX,106HTHE FOLLOWING IS Aw AERQDYN AMICABLY OPTIMIZED ST 
1 AGE CHOSEN FROM THE POSSIBLE FLOWPATH CONFIGURATIONS ABOVE//) 

2224 FQRMAT(//55X, UNSTRESS DATA/SSX , 1 1 H- ///3SX, 4SHALL0WA8LE 

1 AVERAGE DISK STRESS (INPUT) ,F8.0,4H PS I //3SX , 4SHC0MPUTED A 

9 2VERAGE DISK STRESS ,F5.0,4H PSI//35X,45HCGMPUTED R 

200T BLADE STRESS ,F3.0,4H PSI/) 

223S FORMAT (///20X,32HNARN I NG ** NO SOLUTION WHICH SATISFIES THE SPECIF 


1 1 ED DISK STRESS HAS BEEN FOUND ^*///) 

2226 FORMAT ( ///15X , 33hNGTE ... OF 21 SOLUTIONS 08TAINED , IS , 5SH ARE (IS) 
1 WITHIN THE RANGE OF THE SPECIFIED DISK STRESS//) 

2237 FORMAT <20X, 3 1HNOTE PROGRAM PREDICTS STATOR (MEAnLINE) CHOKIN 

1G FOR THE CONDITIONS OF PASS ,12/) 

2223 FORMAT ( 20X , 31HNQTE PROGRAM PREDICTS ROTOR (MEANLINE) ChOKIN 

1G FOR THE CONDITIONS OF PASS ,12/) 

222? FORMAT ( //IX, ' CONTINUITY CANNOT e£ SATISFIED AT THE STATOR EXIT 

IFOR THE ROTOR EXIT SkIRL CONDITIONS SPECIFIED SOLUTION TERMINA 
2TES ' ) 

2240 FORMAT < 59X , 5HSTAGE/59X , SH ///37X , 4SHSTAGE wGRK COEFFICIENT BAS 

1ED ON MEAN RADIUS . , F7 . 3//37X , 4SHSTAGE WORK COEFFICIENT BASED ON H 

2UB RADIUS . . vF7.2//37X,4SHTURBIN£ TOTAl-TOTAl PRESSURE RATIO 

3 »F7. 3/ ) 

2241 F0RMAT(////S5X,11HSTRESS DATA/S5X, 1 1H ///3SX, 45HC0MPUTED 

1 AVERAGE DISK STRESS F3.0,4H PS I //35X , 45HCGMPUTED 

2 BLADE ROOT STRESS ,F3.0,4H PSI/J 

2242 FORMAT (//20X,3SHN0TES ON FLGk ANGLE INPUT //25X,7?h(l) ST 

1AT0R EXIT FLOk ANGLE AND ROTOR EXIT SkIRL CANNOT 5£ SIMULTANEOUSLY 
2 INPUT . /29X, 53HIF BOTH ARE INPUT, ONLY THE ROTOR EXIT SkIRL IS USE 
3D.//25X,32H<2> IF NEITHER ANGLE IS INPUT, STATOR EXIT ANGLE IS CAL 
4CULATED ASSUMING ZERO R0T0R/29X , 10HEXIT Sk IRL//25X , S4H < 3 > IF STATO 
5R EXIT ANGLE IS INPUT, SET ROTOR SkIRL = ZERO (SkIRL wILL BE CAlCU 
6lATED ) / ) 

2243 FORMAT <//20X,96HADDITI0nAu NOTE IF FLGkPATH GEOMETRY IS INPU 

IT, ALLOkABLE DISK STRESS MUST BE INPUT AS ZER0/95X,4h /) 

2244 FORMAT < 20X,72HEACH CONFIGURATION IS CONSTRAINED BY THE FOLLOkIN 
Xii (CONSTANT) PARAME7ERS///35X , 3HRPM/35X , 9HMASS FL0k/35X,44HTUR8INE 
2 INLET TOTAL TEMPERATURE AND PRESSURE/35X, l6rtkORK REGUIREMENT/3SX , 
3 1 5HR0T0R CLEARANCE/3SX ,25HRQTQR EXIT AMIAu MACH NO . /3SX , 1 eriROTGR E 
4XIT SNIRL/3SX ,37HSTAT0R AND ROTOR ZkEIFEL COEFFICIENTS///// ) 

2245 FORMAT ( / ////36X, 13HDISK MATERIAL., 8X, 14H8LADE MATERIAL , SX , 1 7hDISK R 
1IM MATERIAL/24X, 17HDENSITY (LSS/FT3) ,5X, 17HDENSITY (LBS/FT3) ,6X, 

21 7hDENSITY (LSS/FT3) //40X,F5. 1 , 16X, F5. 1 , l?X,F5. 1/ ) 

999 STOP 
END 

SUBROUTINE RCHOKE < VU3 , 0X3 , G , ft ,T1 ,CON,UCHQKE) 

DIMENSION ALE (2) ,ARS(2) , DIFl ( 2 ) , DIF2 ( 2 ) 

ACR°ACRIT ( G , R , Tl ) 

A=lOO. 

S*lOO. 

DEL= 100 , 

K = 0 

20 ALS ( 1 ) =S 
K=K + 1 

IF(K.GE.IOO) GO TO 1400 
VUZ=CON/S+OU3 

T2PP = Tl-M 1 . -(G-l . )/(G+l . )*( <S/AeR)**2>*<2.*(VU2/ACR>/<S/ACR)-l . ) ) 
ACR3=ACRIT(G,R,T2PP) 

ARS ( 1 ) = VU3+SGRT ( ACR3**2- VX3**2 ) 

SS=S+DEL 
ALS(2) =SS 
UU2=C0N/SS+UU3 

T2PP=7l*(l.-(G-l . )/(G+l. )*< <SS/ACR)**2>*<2.*(UU2/ACR>/(SS/ACR>-1 . ) 
1 ) 

ACR3=ACRIT< G ,R ,T2PP) 

ARS ( 2 > o UU3+SGRT ( ACR3**2- 0X3*#2 ) 

DIF1<1)*A8S<AlS<1)-ARS(1>) 

DIFl < 2 > =A0 S ( AlS< 2 ) -ARS ( 2 ) ) 

DIF2(1)=ALS(1)-ARS(1) 

DI F2 ( 2 ) = AlS ( 2 ) -ARS ( 2 ) 

IF(DIF1<1) .LE.5.0) GO TO 1000 
IF(DIFKZ) .LE.5.0) GO TO 1001 
G=DIFZ(1>*DIFZ(2) 

IF(G.LT.O.O) GO TO 200 

_ I F ( D I Ft (.2.) . GT D I F 1 ( 1 ) .) GO TO.. 100 . ... 


200 
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S=S+A*1 .0 
GO TO 20 
DEL»D£L/2. 

A°A/2. 

GO TO 20 
100 D£L=-D£L 
A = -A 

GO TO 20 

1000 UCHQK£ 3 AlS ( 1 > 

GO TO 1200 

1001 UCHOKE°AlS ( 2 > 

GO TO 1200 

1400 UCHOKE 3 AlS< i > 

NRITE(6r 1100) 

1100 F0RMAT(//72H#** ROTOR CHOKE ITERATION LOOP NOT SATISFIED COMPUTA 

1TI0N CONTINUES ***//> 

1200 CONTINUE 
RETURN 
END 

SUBROUTINE DIST < V I N , VOUT » DS » X , VSUC , VPRES s-DUSUC , DUPRES , EXTRM , 
1VRATI0) 

REACT 3 1 . - { V I N/ VOUT ) »*2 
REAC 3 1 .-VIN/VOUT 
G=X-0.S 

AlMAX=0 . 5+EXTRM*REACT 
VMAX=VOUT/< 1 .-DS) 

IF(DS.LT.O.O) GO TO 100 

VCL=>VQUT*( (REAC/2. ) *< SIN < 3 . i4lS?*l2) -1.1 + 1.) 

IF(X.GT.ALMAX)GO TO 10 

VSUC 3 (VI N- VM AX > * ( ( X/ AlMAX ) »*Z) - 2 . * ( VI N-VMAX ) * ( X /ALMAX > + V I N 
DUSUC=>2 . *X/ ( ALflAX**2 ) * ( V I N-VMAX ) -2 . *< VIN-VMAX ) /AlMAX 
GO TO 20 

10 VSUC=VOUT+ ( VQUT-VMAX) / ( ( AlMAX- 1 . ) **2 ) * ( X**2-2 . * AlMAX* < X- 1 , ) -1 . > 
DUSuC 3 ( VOUT- VM AX ) / ( ( ALMAX— 1 . ) **2 > * ( 2 . *X-2 . *AlmAX ) 

20 VPRES 3 VCL- VR AT 1 0* ( VSUC-VCL ) 

DUPR£S=3. 14159*V0UT*( VRATIO+i . ) /2. *REAC*COS (2 . I4l59*li> -VRATIO* 
1DUSUC 
GO TO 40 

100 VCL 3 <VOUT-VIN)*X+VIN 

Al 3 ( VMAX-VOUT ) / ( ALMAX* ( AlmAX-1 . ) >-< VOUT-VIN) /ALMAX 
A2 = ( VOUT-V IN)*( AlMAX+ 1 . ) /AlMAX- ( VMAX-VOUT ) / ( AlMAX* ( AlMAX- 1 . ) > 
VSUC 3 Al*X*X+A2*X+VIN 
DUSUC=2.*A1*X+A2 
VPRES =2.* VCL- VSUC 
DUPR£S=2 .*< VOUT-V I N > -DUSUC 
40 RETURN 
END 

SUBROUT I NE DRAnGE ( G , V IN , VOUT , EXTRM , DSTART , DCHNGE ) 

REACT 3 1 . VIN/VQUT>**2 

XMAX a O. 5+£XTRM*REACT 

VM I N= ( VOUT-V I N ) *XMAX+V I N+0 . 0 1 

VMAX 3 SGRT( (G+l . )/(G-l. ))-0.01 

Dl*l p-VOUT/VMAX 

D2° 1 .-VOUT/VMIN 

DCHNGE 3 ABS ( ( D 1 -D2 ) / 1 0 . > 

DSTART 3 <ABS(D1-D2> >*0.2+D2 

RETURN 

END 

FUNCTION VELIT (WrT»Gl »G2r A lP* AREA ,P,RHO,ACR> 

DIMENSION AlS< 2> rARS(2> *DIF1 (2) ,DIF2<2) 

CONST =H*SGRT ( T/51 S . 7 ) /( COS ( Al P ) *AREA/ 1 44 . * ( P/ 1 4 . 696 ) *RHO*ACR ) 

A=0. 1 
S = 0 . 5 
DEL 3 0.1 
K =0 

20 ALS(1)=S 






K=K + 1 

IF ( K . G£ . 1 00 ) GO TO 1400 

ARS( 1 >=CONST/( ( l .-Gl*(S**2. ) >**G2> 102 

AlS ( 2 ) =S+D£L 

ARS (2) “CONST / ( ( l .-Gl*( (S+D£L)**2. ) )**G2> 

DIFl ( 1 ) =ASS ( ALS < 1 ) -ARS ( 1 ) ) 

DIF? (2) ®A6S ( AlS ( 2 )-ARS ( 2 ) ) 

DIF2<1)=AlS(1)-ARS(1 ) 

DIF2(2)=AlS(2)-ARS(2> 

IF<DIFI ( 1 ) „L£. 0.001 ) GO TO 1000 
IF (DIFl (2) . L£ ,0.001) GO TO 1001 
G®DIF2( l)*DIF2'2) 

IF(G.LT.O.O) GO TO 200 
IF (DIF1(2) .GT.DIFl(l) > GO TO 100 
S®S+A*1 .0 
GO TO 20 
200 D£L®D£L/2. 

A°A/2 . 

GO TO 20 
100 DELa-DEL 
A = -A 

GO TO 20 

1000 VELIT®AlS ( 1 > 

GO TO 1200 

1001 VELIT=AlS<2) 

GO TO 1200 

1400 VELIT=AlS<1) 

WRITE<6» 1100) 

1100 FORMAT ( ///E1H*** SLADE EXIT VELOCITY ITERATION LOOP NOT SATISFIED 
1 COMPUTATION CONTINUES ***///) 

1200' CONTINUE 
RETURN 
END 

SUBROUTINE VNEG ( VI N , VOUT * EXTRM , DS , RATI 0 > 

DIMENSION ALS<2) ,ARS(2) ,DIF2 (2) ,DIF2(2) 

IF(DS.LT.O.O) GO TO 1500 
REACT® 1 .-(VI N / VOUT ) **2 
REAC=1 . -VI N/ VOUT 
ALMAX®0 . 5+EXTRM*R£ACT 
UMAX® VOUT / ( 1 . -DS ) 

CONST® 3 . 1 4 1 59* ( ALMAX**2 ) *VOUT*R£AC/ ( 2 . * ( V I N-VMAX ) ) 

A=0. 1 
S=0.00l 
D£L=0. 1 
K =0 

20 ALS(1)*S 
K®K+1 

IF(K.GE.IOO) GO TO 1400 

ARS < 1 ) =CONST*COS (3.141 59* ( S-0 . 5 > ) +ALMAX 
AlS ( 2 ) =S+D£L 

ARS ( 2 ) =CONST*COS (3.141 59* ( S+DEL-0 . 5 ) ) +AlMAX 
DIFl (1 )®ASS( AlS < 1 )-ARS(I ) ) 

DIFl (2) =ASS ( ALS ( 2)-ARS (2) ) 

DIF2(1)=ALS(1)-ARS(1) 

DIF2(2) =AlS(2)-ARS(2) 

IF(DIFKI).lE.O.OOI) GO TO 1000 
IF(DIF1(2).LE. 0.001) GO TO 1001 
G=DIF2(1 >*DIF2(2) 

IF(G.LT.O.O) GO TO 200 
IF(DIF1(2) . GT .DIFl (1 ) ) GO TO 100 
S=S+A*1 . 0 
GO TO 20 
200 DEL®D£L/2. 

A®A/2. 
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A=-A 

GO TO 20 

1000 XMINoAlS(I) 

GO TO 1200 

1001 XMIN*>AlS<2) 

GO TO 1200 

1400 XMIN=AlS(1) 

NRITE (6,1 100 > 

1100 F0RMAT<///91H*** SLADE PRESSURE SURFACE LOADING ITERATION LOOP NOT 
1 SATISFIED COMPUTATION CONTINUES ***///> 

1200 VCL=VDUT*< (REAC/2. )*CSINC3. 14159* < XMIN-0 . 5) )-l . )+l . ) 

VSUC=<VIN-VMAX)*< ( XMIN/ALMAX ) *-*2 ) -2 . * ( VIN-VMAX) * ( XMIN/AlmAX ) +VIN 

VPRES=2.*VCL-VSUC 

I F < VPRES . GT .0.0) GO TO 1500 

RATIO* ( VCL-0 . 0 1 ) / ( VCL-VPRES ) 

GO TO 1600 
1500 RAT 10* 1.0 
1600 RETURN 
END 

SUBROUTINE AlIFE(BS,TEM,T2ME) 

DIMENSION S< 13) rXi 13) 

C DATA IS FOR IN- 100 MATERIAL 

DATA S/106C0. , 13500. , 17000 ., 21 700 ., 27000 ., 22500 ., 38500 ., 45000 . - 
151500. ,52500. ,66000. ,74000. ,84000. / 

DATA X/52. ,51. ,50. ,49. ,48. ,47. ,46. ,45. ,44. ,43. ,42. ,41 . ,40./ 
IF(BS.L£. 10600. ) GO TO 10 
IF(BS.G£. 84000. ) GO TO 20 
K °1 
KK =2 

40 IF < BS . G£ . S ( K ) . AND .35. l£ . S ( K K ) ) GO- TO 20 

K=K + 1 
KK=KK+1 
GO TO 40 

10 PARAM* ( 52 . —5 1 . )/( 10600.-13500. )*<BS-l3500. )+5l . 

GO TO 50 

20 PARAM=(40.-4i . )/ (84000. -74000. )*<BS-74000. )+4l . 

GO TO 50 

30 PARAM=<X(KK)-X<K> ) / ( S ( KK ) -S < K ) )#<BS-S<K) )+X(K) 

50 G=1000.*PARAM/TEM 

ZFfQ.GE.SO.) TIME=lO.ElO 

IF < G.LT.20. ) TIME* 10 . •>* ( G-20 . ) 

RETURN 

END 

FUNCTION SIMP2 (F , A , DELX , N ) 

C INTEGRATION OF TABULAR FUNCTION, F, BY SIMPSON'S RULE, WHERE 

C U 

C A s lON£R LIMIT, DELX=LENGTH OF SUBINTERVAl, n=NUMBER OF SUBINTERVAl 

C IF N IS LESS THAN 2, SIMP2 IS SET TO 0 

C N-MAY BE ODD OR EVEN 

DIMENSION F ( 51 ) 

23 IF(N-2) 14,17,17 

14 SIMP2=0. 0 
RETURN 

17 IF< <N/2)*2-N) 11,12,14 

12 K*N/2 

ASSIGN 19 TO M 
GO TO 13 

11 K= (N-l )/2 
ASSIGN ie TO M 

13 SUMA=0.0 
DO 15 J=1 ,K 

15 SUMA=SUMA+F(2*J) 

SUMB»0.0 

22 IF(K-l) 14,20,21 

21 DO 16 J=2,K 

16 SUMB=SUMB+F<2*J-l) _ 








20 SIrtP2 3 ':i .3332333*SUr)A+.666e£667*SUr!B+.33332323*F(2*K + l > + .23233332* 

1F< 1) ) *D£LX 

GO TO H,(1?,18) 104 

18 SIMP2 3 SIHPZ+DELX*< . 41&66667*F( Z*K+2 > + .6<se>e>6e>67*F ( Z*K + 1 ) 

1 - . 03232333*F ( 2*K > ) 

19 RETURN 
END 

SUBROUTINE VISCO(T,V> 

Cl=l .264E-5 

C2=0.o00 

TR£F 3 492. 

U 3 Cl*< ( T/TREF ) **C2 > 

RETURN 

END 

SUBROUTINE HM(AJ *FGRM) 

P=l./7. 

TQPSUM b O. 0 
B0TSUM°0 .0 
DO 10 J 3 l*lO 
G 3 J~1 

T0PSU(1 a T0PSUM+ ( < (2.*G+1 . )*(AJ**G> )/< (2.*G+1 . )*P+1. > > 

B0TSUn 3 B0TSUft+( <AJ**G)/< < <2.*G+1 . )*P+1 . )*<2.*<Q+1 . >*P+1 . ) ) ) 

10 CONTINUE 

FORM=TOPSUM/BQTSUM 

RETURN 

END 

SUBROUTINE 5DISK ( RH » RT , CX r RHOD , RHOB , RHOA , S , SOL , A , STRESD r STRESB > 

RH=RH/12. 

RT»RT/12- 

CX»CX/l2. 

A=A/1 44 . 

STRESS 3 (2.51 E-7 ) *RHOS*A* ( S**2 > 

STRESB=1 .2*STRESB 
D£LRA 3 0.5*CX 
B0«0.75*CX 
PHI=0. 1745 
RM=(RH+RT)/2. 

Z=2.*3. 14159*RM*S0L/CX 

RD 3 RH-D£LWA 

BI =80+2 . *RD*TAN ( PHI ) 

SIG1=(3.942E-7)*RH0D*< ( S*RD > **2 > * ( BI /BO+3 . > / ( BI /BO+1 . > 

TBAfi=>0. 1*CX 
ABLD=TBAR*CX 
FB 3 STRESB/1 .2*ASLD*144. 

RCG=RD+D£LRA* ( BO+2 . *CX ) / ( 3 . * ( BO+CX ) > 

0A=3. 14159*RCG*(CX+S0>*D£LRA 

FA=RHQA/ (Z*32. 174)*( (2. *3. 14l59#S/60. )**2>*VA*RCG 
FT =FA+FB 

SIG2 3 Z*FT/ (3 . 14159*RD*80*<BI/B0+1 . ) > 

SIG2 3 SIG2/144. 

stresd 3 Sigi+s:g2 
RH 3 RH*1 2 . 

RT 3 R't + j 2 . 

CX=CX*.vk. 

A=A*144 . 

RETURN 

END 

SUBROUTINE BLOCK (G , T , Z >- BETAI , VI ,BETAG,VQ> 

DIMENSION DIF1(2> ,DIF2(2) r ANGLE (2) ,0CR1 (2) rSS<2> 

Gl=G-l . 

G2=G+1 • 

RAD=3. 14139/180. 

TARGT =B£TAI 
A=RAD 

S 3 8ETAI-S.*RAD 

DEL=RAD ....... 
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K = 0 

20 DO IS 2=1,2 
SS( I )=S 

VACR1=A0S<VI/SIN(S> ) 

• J CRl ( I ) a VACRl 
«FS 1 =G 1 /G2* < VACR 1**2 ) 

vxcri=vacri*cos(e: 

=T/(Z*COS(S) ) 

CM < 1 ,-AFSl *<G2/<Z.*G) ) + <l . -TED* (VXCR 1**2) )/( (1 . -T£T)*VXCRl ) 

D= V I 

E=G*C/G2 

VXCR2=E-SGRT { E**2-l . +G 1 /G2* ( D**2 ) > 

angle < i ) = atan<d/vxcr2) 

VCR2“SGRT < VXCR2**2+D**2) 

IF(I.EG.l) S=S+DEL 
15 CONTINUE 
K = K + 1 

I F ( K . G£ . 1 00 > GOTO 1500 
D I F 1 < 1 ) = ASS < ANGLE i 1 ) -TARGT ) 

DIFl < 2 ) =A8S ( ANGLE < 2 ) -TARGT ) 

DIF2 < 1 ) =ANGLE ( 1 ) -TARGT 
DIF2(2) =AnGLE\2> -TARGT 
IF(DIF1(1) . L£. 0.0000 1 ) GO TO 1000 
IF(DIF1(2).LE. 0.00001) GO TO 1001 
Q=DIF2(1 )*DIF2(2) 

IF(G.lT.O.O) GOTO 200 
IF(DIF1 (2) -GT.DIF1 (1) ) GOTO 100 
GO TO 20 
200 S=S-DEL 
D£L=DEL/2. 

GO TO 20 
100 DEL=-DEL 
S=S+DEL 
GO TO 20 

1000 e£TAO=SS(l) 

VO=VCRl ( 1 ) 

GO TO 1600 

1001 B£TAO=SS ( 2 ) 

VO=VCRl (2) 

GOTO 1600 

1500 BETA0=SS<1> 

VO*VCRl <1 ) 

NRITE < 6 , 1700 ) 

1700 FORMAT ( /// 10X,47H*** BLOCKAGE CALCULATION FAILED TO CONVERGE ***) 
1600 RETURN 
END 

FUNCTION FIT(ZK) 

DIMENSION AK ( 34 ) , AlAm ( 34 ) 

DA i A ALAM / 12. , 1 1 • ,10. , 9 ■ ,3. ,7 .3,7 .6,7 .4,7 .2,7 .052,7 . ,6. 3,6. 6, 6. 

16.2.6. . 5. .4. , j . ,2. ,1. ,0. , — 1 . ,— 2. , — w ■ ,— 4. ,— 5. ,— 6. 7. ,— 3. ,— 9. ,— 1 
2-11 . ,-12./ 

DATA AK /O . 0943 , 0 . 094l ,0.0919,0. 0332 , 0 . 082 1,0.0319,0. 0807 , 0 . 0794 , 

10.0731.0. 0770.0.0767.0.0752.0.0727.0.0721.0.0706.0.0639.0.0599, 

20 . 0497 .0. 0385.0. 0264 . 0 . 0 1 35 . 0 . , -0 . 0 1 40 , -0 . 0234 , -0 . 0429 , -0 . 0575 , 
3-0 . 0720 , -0.0862 , -0 . 0999 ,-0.1130,-0. 1254 ,-0.1 369 ,-0.1 474 ,-0.1 567/ 

IF(2K.G£. 0.0943) GO TO 10 
IF(ZK.L£.-0. 1567) GO TO 20 
K = 1 
KK -2 

40 IF(2K.L£.AK(K) .AND.ZK.GE.AK(KK) ) GO TO 20 
K = K + 1 
KK=KK+1 
GO TO 40 
10 FIT=»l2. 

GO TO 50 

20 FIT=-12.__ _ _ 
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GO TO 50 

FIT°(Ai_AM(KK)-'AlAM(K) >/<AK(KK)-AK(K> >*(2K-AK (X) >+At_AM(K> 

CONTINUE 106 

RETURN 

END 

FUNCTION RECRIT(ZL) 

DIMENSION RE(8),X<8) 

DATA X/8. , 6 . »4. » 2 . » 0 . » -2 • » -4 . » -6 . / 

DATA RE/ 11500. ,9600. ,5600. ,2000. ,645. ,250. , 140. , 100./ 

IF(ZL.GE.a.O) GO TO 10 
IF(ZL.L£.-6.0) GO TO 20 
K a l 
KK=2 

IFtZL.LE.X(K) .A nD.ZL.GE.X(KK) ) GO TO 30 

K = K+1 

KK=KK+1 

GO TO 40 

RECR IT® 11500. 

GO TO SO 

RECRIT® 1 00 . /6 . #ZL+200 . 

GO TO 50 

RECRIT® (RE(KK)-RE(K) ) / { X(KK )-X(K ) )*<ZL-X(K> )+RE(K) 

CONTINUE 

RETURN 

END 

SUBROUT I NE AFLOW ( W , T , P , VX , RHOSTD , ACRSTD , VUH , GAMMA , RH , WCAl , RT ) 
DIMENSION R ( 5 1 ) , FS < 5 1 ) 

VMAX=SGRT < VX**2+VUH**2 ) 

AMAX ■ W»SGRT (T/51S.7)/(P/14. 696 ) / ( RHOSTD* ACRSTD ) / RHGV ( GAMMA , UMAX > 
1*VMAX/VX 

RTMAX=SGRT(AMAX/3.14159+RH**2> 

DELR=(RTMAX~RH>/50. 

KOUNT = l 
RO=RH 
W0 = 0.0 

CONST = 2 . *3 . 1 41 59*VX*RH0STD*ACRSTD* ( P/ 14 . 696 ) /SGRT CT/518 . 7 > 

DO 100 1=1,51 
AJ= 1-1 

R(I>=RO+AI*DELR 

V=SGRT ( VX**2+ ( VUH*RH/R ( I > ) **2 > ‘ 

FS( I > = < ( 1 .-(GAMMA-1 . ) / ( GAMMA+1 . >*<V**2) )**< 1 ./ (GAMMA- 1 . ) ) ) *R ( I > 

CONTINUE 

DO 200 J = 1 ,50 

IF(J.GT.l) GO TO 125 

WCAl®C 0NST* < FS ( 1 ) +FS ( 2 ) ) *D£LR/2 . +wO 

GO TO 150 

WCAL =CONST*S IMP2 ( FS , 0 . 0 , DELR , J ) +W0 
IF(A8S(W-WCAL) .lT. 0.001) GO TO 220 
IF(WCAl.lT.W) GOTO 200 


J J® J-l 

W0=CQNST*SIMP2 < FS , 0. 0, DELR , J J ) +W0 

IF( J J . EG . 1 > WO=CONST* (FS ( 1 ) +FS ( 2 ) ) *D£LR/2. +W0 

DELR=(R( J+l)-R( J))/50. 

RO=R ( J ) 

KOUNT =K0UNT+1 
IF (KOUNT . GE . 5 ) GO TO 260 
GO TO 230 
200 CONTINUE 
260 WR I TE ( 6 » 300 ) WCAL , W 

300 FORMAT ( //5X , 32H *** WARNING *** MASS FLOW ITERATION LOOP DID NOT 

1 CONVERGE ... CALCULATED FLOW = » FS . 4 , 3X , 1 3HACTUAL FLOW =,F8.4/) 

220 RT=R< J + l ) 


RETURN 

END 

FUNCTION RHOV(G,V) 


. RHOV® ( ( 1 , - ( ( G- 1 . > / ( G+l . > )_* < V**2 > > ** U . / < G- 1 . > ) >.* V . . 
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RETURN 
END 

SUBROUTINE SMALSTi ARRAY, LI MIT, MIN) 

DIMENSION ARRAY <200) 

SMALL=ARRAY<1 > 

MIN° 1 

DO 10 1=2, LIMIT 
I F < ARRAY ( I ) . GT . SMAlL > GO TO 10 
SMAlL=ARRAY( I) 

MIN=I 

10 CONTINUE 
RETURN 
END 

FUNCTION AnGIT(H,T,G1 ,G2,VU,AR£A ,P,RH0,ACR,NSLD> 

DIMENSION ALS<2> ,ARS(2> ,DIF1 (2) ,DIF2(2> 
C0NST=W*SGRT(T/5l8.7)/<VU*AREA/144.*<P/i4.696)*RH0*ACR) 

A=0. 0871553 
S=0. 7853976 
AlPMIN=ASIN(VU) 

D£L=0. 0871553 

IF(NBLD.£G. 1 ) DEL=(S-AlPMIN)/10. 

IF ( NBLD . EG , 1 ) A= (S-AlPMIN) /10. 

K=0 

20 ALS(1)=C0S(S)/SIN(S) 

K = K + 1 

IF(K.GE.lOO) GO TO 1400 

ARS ( 1 ) “CONST/ < ( 1 . -Gl* ( < UU/SIN ( S ) )**2 ) ) **G2 ) 

AlS { 2 ) =COS ( S+DEL ) /SIN ( S+DEL ) 

ARS (2 ) “CONST / ( ( 1 . -Gl* ( ( UU/SIN (S+DEL ) >**2) )**G2) 

DIF1 < 1 ) =A8S < AlS < 1 ) -ARS < 1 ) ) 

DIFl i 2 ) = ABS ( AlS ( 2 > -ARS < 2 ) > 

DIF2<1)=AlS<1)-ARS(1) 

D I F2 ( 2 ) = ALS < 2 ) -ARS < 2 ) 

IF(DIF1(1 ) .lE. 0.002) GO TO 1000 
IF (DIFl (2) .LE. 0.002) GO TO 1001 
Q=DIF2( 1 )#DIF2<2) 

IF(G.LT.O.O) GO TO 200 
IF(DIF1(2) .GT.DIFl(l) ) GO TO 100 
S=S+A*1 .0 
GO TO 20 
200 DEL=D£L/2. 

A=A/2 . 

GO TO 20 
100 D£L=-DEL 
A=-A 
GO TO 20 

1000 BNGIT=AlS(1> 

ANGIT=ATAN( 1 . /BNGIT) 

GO TG 1200 

1001 BNGIT=AlS (2) 

AnGIT=ATAn( 1 . /BNGIT) 

GO TO 1200 

1400 BNGIT=AlS(1> 

ANGIT=ATAN( 1 ./BNGIT) 

WRITE(6, 1500) 

1500 FORMAT (///15X.77H*** BLADE EXIT AnGLE ITERATION LOOP NOT SATISFIED 
1 COMPUTATION CONTINUES ***///) 

1200 CONTINUE 
RETURN 
END 

FUNCTION ACRIT (G,R,T) 

ACR I T=SGRT ( 2 . 0*G/ < G+l . ) *32 . 1 74*R*T ) 

RETURN 

END 

SUBROUTINE FD I SK ( GAMMA , VEX I T , GAMM 1 , GAMM2 , GAMM3 , NPASS , K PASS , RGAS , 

1 SPEED , P 1 P , T 1 P , RNORK , NFLGN , ALP 1 , RM 1 , RM2 , RM3 , Hi , H2 , H3 , T2P , U2 , UACR2 , 
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2U3 , V 1 , VU 1 » ETA t P2P r~3P r P2P r PRTURB , AlP 2 , V2 r VU2 , U'*RK2 , VHRK3 , UACR3 , 
3VU2 , TPTP2 , T2PP , PDPP2 , 732PP , 73PP * i*U2, HUPP2 » SE7A2 r WPP2 , WU3 , WUPP3 , V2 , 
4 AlP 3 f UX3 1 HXPP3 r WPP3 , SET A3 r P32PP , PLOS 1 , PL0S2 , RH3 , RT3 r RTS » ACRl r ADP3 •- 
5H1 A» rt2 A , m 3A r VX2 , HFACTM , WFACTH , UT3 , STGACC , RSEST » 72PPH , ACR3 ) 
PASS=KPASS 
T2P°TlP 
TSTD a 518. 7 

* ASTD a ACRTT (GAMMA rRGASrTSTD) 

i RHSTD a 2ll6.22/<RGAS*TSTD) 

CP= ( 1 . /GAP1M2 ) *RGAS/778 . 26 
1 ACRl=ACRIT (GAMMA»RGAS»71P) 

‘ T3P a Tl P-RHORK /CP 

: ACR3»ACRIT (GAMMA ,RGAS »73P) 

* UMAX**ACR3 

i' VX3 a V£XIT 

! C0N«25039 . 737*RHQRK 

’ UX3D=VX3*ACR3 

( U3 a UX3/COS(ALP3) 

i MU2»VX3 *TAn<AlP3) 

\ VWRK3=UU3*ACR3 

, IF(KPASS.GT.O) GO TO 10 

IF(RBEST.EG.O.O) CAlL RCHOKEi UHRK3 r UX3D, GAMMA ,RGAS ,Tl P ,CO(M # UCHOKE ) 
| IF ( UCHOKE. LT.ACR3) UMAX=UCHQK£ 

IF(AlP1 .LT. 0.349) GO TO 100 

| if( alpi . eq . Alps) go to 100 

| IF i ALP3 .GT .0.0) GO TO 100 

UFLA7 a COlV/ < UX3D* < TAn ( Al_Pl +0 . 140)-7 AN<AlP3> ) ) 

* IF <UFLAT.LT. UMAX) UMAX=UFLAT 

100 IF<R8EST.£G.0.0)G0 TO 20 

RM2 a RBEST 

U3=RM3*3. l4159*SPEED/30. 

! GO TO 15 

20' UMIN a <VWRK3+SGRT(VU*RK3**2+2.*C0N) ) /2. + 10. 

30 RMIN=UMIM/ ( 3 . l4l59*SP£ED/30 . ) 

D£L a ( UMAX-UMIN) /20. 

U3=UMIW 
RM3 a RMlN 
GO TO 15 

10 U3“UMIN+PASS*DEL 

S RM3 a U3/ <3. l4l59*SPE£D/30. ) 

| 15 SWPM=<U3**2)/<25036.62*RW0RK) 

HFACTM=1 ./SHPM 

IF<AlP2.EG.O.O) STGACC=1 .-WFACTM/2. 

IF ( NPASS. EG . 1 ) ETA= (0. 92-&SHPM ) / (SMPM+O . 0227 ) 
IF(NPASS.EG.l)P2P=0.98*PlP 

P3P a PlP*< < 1 .-( 1 . -73P/T1 P ) /£7A ) ** ( 1 . /GAmM2 ) ) 

PRTURS a Pl P/P3P 
UACR3°U3/ACR3 

AREA3 a HFLQW#SGRT ( T3P/TSTD ) / < P3P/ 1 4 . 696 ) / ( RHSTD*ASTD ) / f RHGV C GAMMA , 
1V3)*C0 S(AlP3) ) 

H3=AR£A3/(2.*3. 14159*RM3) 

H3A»H3 

RH3=RM3-H3/2. 

I F < RBEST . GT .0.0) GO TO 25 
! I F < K PASS . GT . 0 ) GO TO 25 

IF<RH3.G£. 0.083) GO TO 25 
UMIN a UMIN+lO. 

GO TO 30 

; 25 RT3=>RH3+H3 

VWRK2=239lll .88#RH0RK/(SPEED*RM3)+VHRK3 

UU2 = VWRK2 / ACR 1 

UX2D a UX3D 

VX2=MX2D/ACR1 

) V2 a SGRT(VU2**2+UX2**2) 

ALP2 a AT AM ( VU2/ VX2 ) 

U2=U3 _ _ . _ . _ 
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0ACR2=U2/ACRl 

> RH2 = Rfl3 
RH2=RH3 

HFACTH=HFACTM* ( ( RN3/RH3 ) **2 ) 

UT3=U3*RT3/Rtt3 

VUH2»RM2*VU2/RH2 

CALL AFLOH ( NFLOH , T2P » F2P » UX2 , RHSTD , ASTD , UUh2 , GAfinA > Rh2 , wCAl , RT2 ) 
t. H2A=RT2-RH2 

H2=WFLOH*SGRT < T 1 P/TSTD >/< < RHOU < GAWMA , V2 ) *VX2/ V2 ) *2 . *3 . 1 4 1 59*RM2* 
I P2P/ 1 4 . 69o*RHSTD*ASTD ) 

RTS=RT2 

AREA 1=3.141 59* ( RT2**2-RH2*#2 ) # 1 44 . 

SLPl =ASS i ALPl ) 

k V 1 a VEL I T ( HFLOH » T 1 P » GAfJH 1 , GAMt*J3r SLPl ,AREAl , PIP, RHSTD, ASTD) 

VU1=V1*SIN(ALP1> 

Hi A=h2A 
RM1 =RM2 

Hi =wFLQW*SGRT ( T1 P/TSTD ) / ( (RHCV ( GAMMA »Ul ) *COS ( A lP 1 ) ) *2 . *3 . 14159* 
1RM1*P1P/14.696*RHSTD*ASTD> 
j TPTP2 3 1 . -GAClMl * ( UACR2**2 ) * ( 2 . *VU2/UACR2-1 . ) 

> T2PP=TPTP2*T2P 
UACRH2=RH2/RM2*UACR2 

1 TPTPH2=1 .-GAnm*<UACRH2**2)*<2.*VUH2/UACRH2-l . > 

j T2PPH=TPTPH2*T2P 

| PDPP2 = TPTP2**< 1 . /GAWJ2> 

! T32PP=1 . ~G AMM 1 * < UACR2**2 ) * < 1 ./TPTP2)*( 1 <RM3/RM2)**2) > 

I T3PP»T32PP*T2PP 

1 ADP2=ACR I T ( GAP?MA i RGAS r T2PP > 

J ADP3=ADP2*SGRT(T32PP> 

l*U2 3 VURK2-U2 
HUPP2=WU2/ADP2 

BETA2=ASIN< i VU2-UACR2 ) / ( SGRT < V2**2+UACR2**2-2 . *VU2*UACR2 ) ) ) 

) HPP2 3 UUPP2/SIN<8£TA2> 

WU3=VHRK3-U3 

HUPP3=WU3/ADP3 

WXPP3*VX3*ACR3/ADP3 

HPP3=SQR7(HXPP3#*2+WUPP3**2) 

W3FPS=WPP3#ADP3 
I H2FPS 3 WPP2*ADP2 

| V2FPS=V2*ACRl 

j V3FPS=V3*ACR3 

! IF ( ALPS . ME . 0 . 0 ) STGACC= ( H3FPS**2-H2FPS**2 ) / i V2FPS**2-U2FPS*#2+ 

1H3FPS**2-U2FPS**2 ) 

BE7A3=ASIN < HUPP3/HPP3 ) 

v IF( NPASS . EG ■ 1 ) P32PP=T32PP** ( 1 . /GAMM2 ) 

f IF ( NPASS. EG . 1 ) PLOSl =0 . 9B 

IF<NPASS .EG . 1 ) PL0S2=P32PP 
RETURN 
END 
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